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A Complete Self-Contained Audio Measurement System, by James D. Foote This new 
audio analyzer has everything needed lor audio measurements — source, filters, detectors, 
voltmeter, and counter — all under microprocessor control. 

Audio Analyzer Applications The major areas are general audio testing, trans- 
ceiver testing, and automatic systems. 

Making the Most of a Microprocessor-Based Instrument Controller, by Corydon J. Boyan 

In an audio analyzer, microprocessor control means automatic operation, "guaranteed" 
accurate measurements, and extra features. 



Design for a Low-Distortion, Fast-Settling Source, by George D. Pontis It's based on a 
state-variable filter with refinements. 

Floating a Source Output, by George D. Pontis The floating output lets the user eliminate 
ground-loop errors, sum signals, and add dc offsets. 

A Digitally Tuned Notch Filter, by Chung Y. Lau It eliminates the fundamental frequency 
component of the incoming signal for distortion and noise measurements. 

A Custom LSI Approach to a Personal Computer, by Todd R. Lynch Nine HP-produced 
large-scale integrated circuits make the HP-85 possible. 

Handheld Calculator Evaluates Integrals, by William M. Kahan Now you can carry in your 
pocket a powerful numerical integrator like those available on large computers. 



In this Issue: 

Audio frequencies are frequencies within the range of human hearing, roughly 20 Hz to 20 
I kHz. Frequencies on either side of this range are often loosely classified as "audio,'' too. 
I Thus, the frequency range of the instrument shown on the cover, Model 8903A Audio 
Analyzer (page 3). is 20 Hz to 100 kHz. The 8903A is used for testing— among other 
^J-Jj^J ~— things — many of the electronic devices that speak to us or play music, such as CB radios and 
^^•^^^^^^ high-fidelity systems. In our cover photo it's shown plotting the frequency response of the 
stereo amplifier on the left at different signal levels. 
There are, of course, other ways of making the measurements in the 8903A's repertoire. 
What makes the 8903A better? First, it's a complete system that includes a low-distortion signal source, a 
counter, a voltmeter, and various filters and detectors. Second, all of this is under microprocessor control, 
automatically stepping through complicated sequences of measurements and computations. Third, it's ex- 
tremely accurate; for example, it can measure total harmonic distortion (THD) down to 0.003% under normal 
conditions. Fourth, the 8903A has recorder outputs that make plotting results about as easy as it can be. 

If you're not so fortunate as to have studied integral calculus in school, the integral of a function probably isn't 
the familiar and highly useful concept that it is to scientists and engineers. One way to think of an integral is this: 
draw the graph of the function as a meandering line on a piece of graph paper. Then the integral of the function is 
the area bounded by 1) the graph of the function, 2) the horizontal axis of the graph paper, and 3) two vertical 
lines called the upper and lower limits of integration. Sometimes the integral of the function can be expressed 
neatly in mathematical terms, but more often than not it can't. So various methods have been devised for 
estimating integrals using computers. Because most of these numerical integration programs run on very, very 
large computers, it seems like a small miracle that you can now carry a numerical integrator — a very good 
one — around in your pocket. Beginning on page 23. its designer tells us about its capabilities and limitations. 

On page 18 is an article about the nine special integrated circuits that make the HP-85 Computer possible. 
This set of custom integrated circuit chips minimizes the cost of the electronics, eases the problem of cooling the 
computer, makes the small package possible, and provides features that couldn't have been included otherwise. 
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A Complete Self-Contained Audio 
Measurement System 

This automatic, autoranging audio analyzer has the signal 
source, distortion analyzer, and counter to make the 
measurements most often needed in audio-frequency 
testing. 

by James D. Foote 



HEWLETT-PACKARD'S NEW MODEL 8903A Audio 
Analyzer (Fig. 1) is a complete audio measurement 
system for quick and accurate characterization 
of systems and signals in the frequency range 20 Hz to 100 
kHz. The starling point for the 8903A is the classical distor- 
tion analyzer. Added to this are microprocessor control, a 
reciprocal frequency counter, rms detectors, and a pro- 
grammable audio source. These provide accurate mea- 
surement of ac level, distortion. S1NAD," signal-to-noise 
ratio, and dr. level. The audio source and the measurement 
circuits can work independently or together. The source is 
programmable in frequency and level and has very low 
distortion. The measurement circuits can monitor this in- 
ternal source or any other independent input waveform. 
Together the source and measurement input can be used for 
swept response measurements. 

'SiNAD = 'alio ol signal plus noise plus distortion 10 noise plus distortion in a receiver s 
output 



All measurements are available at the push of a button. 
No knob adjustment or operator interaction is necessary. 
One simply applies the signal and selects the measurement 
mode. All control and processing are handled by the inter- 
nal microprocessor. The microprocessor monitors the input 
signal and makes internal gain and frequency adjustments 
as required. 

In automatic measurement systems, the 8903 A is capable 
of rapid and straightforward remote control. Analyzer op- 
erations can be controlled and all measurements can be 
transferred via the Hewlett-Packard Interface Bus (HP-IB), 
Hewlett-Packard's implementation of IEEE Standard 488- 
1978. On the bench, the 8903A allows rapid and accurate 
circuit characterization when many repetitive measure- 
ments are necessary. 

Major application areas for the 8903 A Audio Analyzer are 
general audio testing, transceiver testing, and automatic 
systems. In general audio testing, the 8903A measures the 




c 



33 32 " J ~ 



J Tl 




Fig. 1. Model 8903 A Audio 
Analyzer makes the accurate 
measurements needed to charac 
tenze systems and signals in Us 
frequency range ol 20 Hz to 100 
kHz It has applications m general 
audio testing, transceiver testing, 
and automatic systems Micro- 
processor control makes it auto- 
matic and easy to use. 
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Fig. 2. The 8903A Audio Analyzer is basically a distortion analyzer, with a tunable notch lilter to 
remove the tundamental frequency component of the signal and a detector to measure what 
remains, which consists of noise and distortion Added to this are a microprocessor-based 
controller, a reciprocal frequency counter, a programmable audio source, and other resources, 
forming a complete audio measurement system. 



frequency response and distortion characteristics of filters, 
high-quality amplifiers, audio integrated circuits, and other 
devices. The frequency of the internal source can be swept 
while making measurements in all modes. The analyzer 
provides recorder outputs and scaling for easy generation of 
plots using an X-Y recorder. 

For transceiver applications the most common receiver 
measurements are SINAD for FM receivers and signal-to- 
noise ratio for AM receivers. A psophometric filter is in- 
cluded for making measurements to CEPT standards. 
Common transmitter measurements such as audio distor- 
tion and squelch tones are made using the 8903A with its 
companion instrument, the 8901 A Modulation Analyzer. 1 
In automatic systems, the 8903A provides many frequently 
needed audio functions, doing the work of an audio syn- 
thesizer, digital multimeter, frequency counter, and tun- 
able notch lilter. More details on specific applications are 
presented on page fi. 

Control Philosophy 

Front-panel control of the audio analyzer is simple, yet 
powerful. Most functions can be used and understood with 
very little training. The casual user can select amplitude, 
frequency, measurement mode, and filtering simply by 
reading the labels on Ihe controls. More details are available 
on the instrument's pull-out card. 

A great deal of measurement sophistication is built into 
the 8903A software. Measurement routines are structured 
to optimize measurement speed and accuracy. As a rule, 



measurements triggered from the bus or initiated from the 
keyboard are accurate from the first reading. The operator 
needn't wait for successive measurements to verify that the 
reading has stabilized. The software algorithm monitors 
key voltages in the audio chain and waits until they 
stabilize before taking data. Not only does the software 
perform these functions much more rapidly than the 
operator, but it can also ensure the optimal convergence of 
the measurement with a repeatable, well defined technique. 
Distortion. SINAD. and signal-to-noise ratio, in particular, 
are examples of measurements that in the past required a 
significant amount of settling time and operator interac- 
tion. A classical distortion analyzer requires repeated ad- 
justments to achieve an accurate distortion reading. More 
recent analyzers have offered semiautomatic tuning and 
leveling, but response time is often long, and operator in- 
teraction is required if the frequency, amplitude or relative 
distortion of the signal changes significantly. With the 
8903 A. some delay still exists, but the delay involved is 
minimized by careful circuit design and microprocessor 
control. 

Special functions extend user control of the instrument 
beyond thai normally available from the front panel. These 
functions are intended for the user who knows the instru- 
ment and the service technician who needs arbitrary con- 
trol of the instrument functions. Automatic tuning and 
ranging, overvoltage protection, and error messages protect 
the user from invalid measurements during normal opera- 
tion. When special functions are used, some of these 
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safeguards are removed, depending on the special function 
selected, and thus there is a degree of risk that the mea- 
surement may be invalid. However, there is no risk of dam- 
age to the instrument. 

To enter a special function, the user enters the special 
function code (usually a prefix, decimal, and suffix) then 
presses the SPCL key. The special function code appears on 
the display as it is being entered. If a mistake is made during 
entry of the special function code, the user can press the 
CLEAR key and start over. When a special function is en- 
tered, the light in the SPCL key goes on if it is not already on. 
The readout on the display depends on the special function 
entered. It may be a measured quantity, an instrument set- 
ting, or a special code. In some cases the display is unal- 
tered. Special functions can be entered from the HP-1H by 
issuing the special function code followed by the code SP. 

Floating Input and Output 

To eliminate troublesome ground loops, both the source 
output and the measurement input of the 890:iA are float- 
ing. This is helpful in low-distortion or low-level ac mea- 
surements when it is necessary to reject potential differ- 
ences between the chassis of the 8903A and the device 
under test. The 8903A also has EMI (electromagnetic inter- 
ference) protection built into the source output and mea- 
surement input lines so that it can work in the presence of 
high EMI. All of the analog circuitry is shielded by an 
internal EMI-tight box. The output and input lines extend- 
ing from this box to the front panel are shielded and termi- 
nated in BNC connectors. For user convenience, the BNC 
connectors are spaced so that BNC-to-banana adapters can 
be attached. Thus a banana or twisted-wire connection can be 
made to the instrument when EMI shielding is not critical. 

Analyzer Architecture 

The 890.IA Audio Analyzer combines three instruments 
into one: a low-distortion audio source, a general-purpose 
voltmeter with a tunable notch filter at the input, and a 
frequency counter. Measurements are managed by the 
microprocessor-based controller. This combination can 
make most common measurements on audio circuits au- 
tomatically. To add to its versatility, the analyzer also has 
selectable input filters, logarithmic frequency sweep, X and 
Y outputs for plotting measurement results versus fre- 
quency, and HP-IH programmability. Fig. 2 is a simplified 
block diagram. 

The amplitude measurement path flows from the INPUT 
jacks (HIGH and LOW) to the MONITOR output on the rear 
panel, and includes the input and output rms detectors, the 
dc voltmeter (the voltage-to-time converter and counter), 
anil the SINAD meter circuitry. Measurements are made on 
the difference between the signals at the HIGH jack and the 
LOW jack. Differential levels can be as high as 31)0 V. Signals 
that are common to both the HIGH and LOW jacks are bal- 
anced out. Signals applied to the LOW jack must not ex- 
ceed 4 V. 

The input signal isac coupled for all measurement modes 
except dc level. The signal is scaled by the input attenuator 
to a level of 3V or less. To protect the active circuits, the 
nvervoltage protection circuit quickly disconnects the 
input amplifier if its input exceeds 15V, 



The differential signal is converted to a single-ended 
signal (referenced to ground) and amplified. The signal is 
further amplified by a programmable gain amplifier, which 
is ac coupled. The gain of this amplifier and the 
differential-to-single-ended amplifier are set to keep the 
signal level at the input rms detector between 1.7 and 3V 
rms to optimize its effectiveness and accuracy. 

The output from the first programmable gain amplifier is 
converted to dc by the input rms detector and measured by 
the dc voltmeter. The output of the detector is used to set the 
gain of the input circuits and becomes the numerator of the 
SINAD measurement and the denominator of the distortion 
measurement. The gain of the input path is determined by 
measuring the dc level. The input rms detector also 
monitors the ac component (if there is one) and lowers the 
gain of the input path if theacsignal will overload the input 
amplifier. At this point either the 400-Hz high-pass filter or 
the psophometric filter can be inserted into the signal path. 
The 400-Hz high-pass filter is often used to suppress line 
hum or the low-frequency squelch tone used in some 
mobile receivers. The psophometric filter has a bandpass 
frequency response that simulates the "average" response of 
human hearing. It is often used to condition a receiver audio 
output when determining the receiver's input sensitivity. 

During SINAD. distortion, or distortion level measure- 
ments, the fundamental of the signal is removed by the 
notch filter. The output from the filter is the distortion and 
noise of the signal. In the ac level and signal-to-noise 
modes the notch filter is bypassed. After amplifying and 
low-pass filtering, the output from the notch filter is con- 
verted to dc by the output rms detector and measured by 
the dc voltmeter. 

During distortion or distortion level measurements, the 
notch filter is tuned to the frequency counted at its input. 
Coarse tuning is done by the controller, and internal analog 
circuitry fine tunes and balances the notch filter. During 
SINAD measurements, the controller coarse tunes the notch 
In the source frequency. Thus a SINAD measurement is 
normally made with the internal source .is the stimulus; this 
permits measurements in the presence of large amounts of 
noise (where the controller would he unable to determine 
the input frequency). If an external source is used in the 
SINAD measurement mode, the source frequency must 
be within 5% of the frequency of the internal source. 

The two programmable gain amplifiers following the 
notch filter amplify the low-level noise and distortion sig- 
nals from the notch filter. The overall gain of the two 
amplifiers is normally set to maintain a signal level of 0.25 
to 3V at the output detector and monitor. The 30-kHz and 
80-kIIz low-pass filters are selected from the keyboard. 
With no low-pass filtering, the bandwidth of the measure- 
ment system is 750 kHz. The filters are most often used to 
remove the high-frequency noise components in low- 
frequency distortion and signal-to-noise measurements. 
The output Irom the second programmable gain amplifier 
drives the rear-panel MONITOR output jack. Taking advan- 
tage of the increased amplification available at this point, 
the counter monitors this output in ac level and signal-to- 
noise modes. 

The output rms detector is read by the dc voltmeter in the 
ac level, SINAD (the denominator), distortion (the 
numerator), distortion level, and signal-to-noise measure- 
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Audio Analyzer Applications 



The 8903A Audio Analyzer's measurement capabilities reach far 
beyond conventional distortion analyzers. Much of this performance 
results from microprocessor control and HP-IB programmability 
Numerous hardware features such as a fast counter, both analog 
meter and digital display, and switchable detector filtering allow the 
user to make unusual or special measurements with convenience 
and little auxiliary apparatus 
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men! modes. It is also used to set the gain of Ihe two pro- 
grammable gain amplifiers. The detector Can be configured 
internally to respond to the average absolute value of the 
signal instead of the true rms value. This option is provided 
because some measurement Specifications for detection of 
distortion and noise specify the use of an average respond- 
ing detector. Average responding detectors do not give an 
accurate indication of signal power unless the signal wave- 
form is known. (If the waveform is Gaussian noise the read- 
ing will be approximately 1 (IB less than the true rms value.) 

In the SINAD mode the outputs from the input and output 
rms detectors are converted to logarithms, subtracted, and 
converted to a current by the SIN'AD meter amplifier to drive 
the SINAD panel meter. Since SINAI) measurements are 
often made under very noisy conditions, the panel meter 
makes it easier to average the reading and to discern trends. 

The voltage-to-time converter converts the dc inputs into a 
time interval, which is measured by the counter. 

The 8903A uses a reciprocal counter. To measure fre- 
quency, it counts the period of one or more cycles of the 
signal at its input, Then the controller divides the number of 
cycles by the accumulated count. The reference for the 
counter is the 2-MHz time base, which also is the clock for the 
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Consider the 8903A used at a test and calibration station in the 
manufacture of audio power amplifiers A typical sequence of events 
might include an output offset null, Irequency response check, distor- 
tion test, and noise measurement The 8903A can perlorm all these 
measurements quickly with a single test setup It an X-Y plotter is 
connected to the rear-panel outputs, the results of swept frequency 
measurements can be recorded on standard log/log or log/tin 
graph paper 



Fig. 1. Swept distortion and Ire- 
40 kHz 50 kHz 60 kHz quency response of a two-pole ac- 

tive Utter, measured by the 8903A 
Audio Analyzer. 



controller. The counter has four inputs and three modes of 
operation: 

1. Voltage measurement. The time interval from the 
voltage-to-time converter is counted. The accumulated 
count is proportional to the dc voltage. For direct mea- 
surements |ac level and distortion level), the count is 
processed directly by the controller and the result is dis- 
played. For ratio measurements (SINAD, distortion, and 
signal-to-noise], the counts of two successive measure- 
ments are processed and displayed. For SINAD and dis- 
tortion, the controller computes the ratio of the outputs of 
the input and output rms detectors. For signal-to-noise 
measurements, the output of the output rms detector is 
measured with the oscillator on and off and the ratio of the 
two measurements is computed. 

2. Input frequency measurement. The signal from the last 
programmable gain amplifier or the high-pass bandpass 
filters is conditioned by the counter input Schmitt trig- 
ger to make it compatible with the counter's input. The 
period of the signal is then counted, the count is pro- 
cessed by the controller, and the frequency is displayed. 

3. Source frequency measurement. The counter measures 
the frequency of the oscillator during tuning and when 
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Fig. 2. Test setup lor screening operational amplifiers 

The swept frequency measurement with plot capability finds many 
applications in the laboratory Fig 1 shows the swept distortion and 
frequency response of a two-pole active filter The upper curve shows 
the magnitude response, while the lower shows distortion. Notice that 
the analyzer input magnitude covers a 30-dB dynamic range. During 
the sweep the 8903A is automatically setting input gam before per- 
forming the distortion measurement 

If an HP-IB controller is available, there are many more applications 
for the 8903A Fig. 2 shows a simple test setup for screening opera- 
tional amplifiers With no other instruments in the system, a 
computer-controlled 8903A can rapidly and accurately measure 
input offset voltage, input noise voltage, and distortion It can also be 
used to measure the gain-bandwidth product of the op-amp, pro- 
vided it is not greater than 30 MHz The controller test program can be 
written to provide either a go/ no-go output or a listing of measurement 
results 

Many of the special measurement modes are available through the 
use of special lunctions. For example, the 8903A can be used as a 
test amplifier with gain settable Irom -24 dB to +94 dB. Signal 
filtering can be added by selecting the appropriate front-panel con- 
trols, or the special functions can be used to put the instrument into a 
notch or bandpass filter mode Of course, signal Irequency and 
amplitude will be measured and displayed if the operating loads are 
chosen correctly. This mode of operation makes it possible to count 
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Fig. 3. Transmitter test setup using the 8901 A Modulation 
Analyzer and the 8903A Audio Analyzer. 

the frequency of a signal whose amplitude is m the low tens of 
microvolts. 

Transceiver Testing 

Much effort has gone into the design of the 8903A to facilitate the 
audio measurements required for automatic programmed, 
transceiver testing. For example the worst-case source frequency 
error of 0.3% allows the 8903A source to replace a synthesizer for 
squelch-tone generation in addition, a binary programming mode is 
available through the HP-IB that causes the 8903A to generate a tone 
burst sequence that can be used to unlock a coded receiver A 
related function allows the 8903A to measure burst tones generated 
by an external source, such as a transmitter under test The packed 
four-byte output allows the 8903A to output frequency measurements 
as often as every eight milliseconds 

When the 8903A is used in coniunction with the 8901 A Modulation 
Analyzer, almost all transmitter tests can be automated. Fig. 3 shows 
the block diagram. With the source turned oft and the transmitter 
keyed, the squelch-tone frequency can be counted Then, the 
400-Hz high-pass filter can be switched in to eliminate the squelcn 
tone from the remaining measurements With the source output turned 
on. the 8903A can easily be programmed to make the necessary 
measurements to determine distortion and microphone sensitivity 



verifying that the oscillator frequency is within toler- 
ance. This frequency is normally not displayed. 

The source covers the frequency range of 20 Hz In 100 
kHz. It is tuned by the controller to the frequency entered 
from the keyboard, using a tune-and-count routine. (Note 
that the frequency is not obtained by frequency synthesis.) 
The switch following the oscillator is closed except in the 
signal-to-noise measurement mode or when an amplitude 
of OV is entered from the keyboard. 

The source output amplifier and output attenuators pro- 
vide 77.5 dB of attenuation in 2.5-dB steps, and 2.5 dB of 
attenuation in 250 steps. This gives an open-circuit output 
from 0.6 mV to 6V. The floating output amplifier converts 
the ground-referenced input to a floating output. Hither 
output. HIGH or LOW, can be floated up to ten volts peak. 

The entire operation of the instrument is under control of 
the microprocessor-based controller, which sets up the in- 
strument at turn-on. interprets keyboard entries, executes 
changes in the mode of operation, continually monitors 
instrument operation, sends measurement results and er- 
rors to the front-panel displays, and interfaces with the 
lll'-IB. lis computing capability is also used to simplify 
circuit operation. For example, it forms the last stage of the 



counter, converts measurement results into ratios ( in % or 
dB|.and soon. It also executes routines for servicing the rest 
of the instrument as well as itself. 

Input Circuits 

Numerous design constraints were imposed on the input 
attenuator protection'amplifier circuitry. An input imped- 
ance of 100 kSl is necessary to prevent the input circuitry 
from unnecessarily loading the device under test. On the 
other hand, maintaining a good signal-to-noise and distor- 
tion ratio, good frequency response to 100 kHz, automatic 
operation and reliable performance with input signals as 
large as 300V is very challenging. Consider the 300V con- 
straint. It is possible for a 300V signal to appear suddenly at 
the input while the instrument is measuring a 50-mV input 
level. Not only must no damage occur, but also the overload 
recovery must be quick, and the input protection circuit 
must not be allowed to degrade the input noise floor when 
measuring 50 mV. The 300V signal may also have spikes or 
transients that far exceed 300V, or the user may inadver- 
tently apply a larger signal. In all these cases the circuit 
must recover without causing a safety hazard to the user, 
destroying internal components, or even blowing a fuse. 
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Making the Most of a Microprocessor-Based 
Instrument Controller 



by Corydon J. Boyan 



In the 8903A Audio Anaiy£er. most of the tasks that lead to the 
display ot measurements are coordinated by a microprocessor The 
processor (an F8 with 18K bytes of ROM and 192 bytes of RAM) 
counts and tunes the internal audio source, sets the input amplifier's 
gam. tunes the notch filter, sets the output amplifier's gain, and 
measures the voltages that will be used to generate each reading 
This means that the 8903A can. among other things, automatically 
take distortion readings without calling upon the user to turn several 
knobs when seeking a null This ability alone is ample |ustitication for 
basing the instrument controller on a microprocessor, but the 8903A 
goes far beyond this in applying the power of its controller 

"Guaranteed" Accurate Measurements 

The 8903A is HP-IB programmable, and this brought some impor- 
tant factors into consideration during the design. The performance of 
the internal source when settling from one frequency or level to 
another, for example, is a key tactor in assuring the validity of the first 
measurement taken after changing the frequency or level Similarly, 
the settling performance of the input and output amplifiers and the 
notch filter under various adverse signal conditions becomes very 
important when the instrument is trying to deliver an accurate first 
measurement after a change in operating parameters (e.g., after 
tuning to a new frequency) Add to this the desire to make measure- 
ments as rapidly as possible and you have a very interesting problem 
for a microprocessor to solve. 

For example, every time the internal source is tuned, the processor 
spends about 1 70 milliseconds tuning to and verifying the frequency 
This translates to over 86,000 operations. Of equal complexity is the 
|ob of setting the correct gam for the output amplifier and allowing the 
circuits to settle before the output amplifier voltage is read, with the 
object of never giving the user an invalid reading while at the same 
time delivering the reading in as short a time as possible To ac- 
complish this, the routine that controls the output amplifier makes use 
of such techniques as slope checking and frequency dependent 
delays to ensure rapid, valid readings These techniques make pos- 
sible accurate readings in half the time it might otherwise take. 
Referring to the flow chart, Fig 1 , note that the key to this routine is the 
technique of measuring the rate of change of the signal on the output 
amplifier (the slope) before attempting to check the signal level to 
determine if the output amplifier gam is properly set. This is done 
because it is quite common, after the signal level has suddenly 
changed (for example, when the notch filter is suddenly switched in 
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Fig. 1. 8903A leveling algorithm. The microprocessor checks 
the rale ol change ol Ihe input signal as well as its level to make 
certain that the signal has settled before the output amplifier 
gam is set. 

when going from ac level fo distortion mode), to have a rapidly falling 
output amplifier voltage pass thiough the acceptable region, and 
thus have the gam appear to be properly set. when in fact the signal is 
on its way to a level that will require more gam for the proper reading 
To keep the controller from being fooled by this phenomenon, we take 
two voltage readings in rapid succession, and from them calculate 
the rate of change of Ihe signal. The time over which we measure this 
slope is longer at lower frequencies, and is m fact looked up in a table 
m ROM, based on the frequency at which we're operating If the rate 
of change is too fast during this period, we delay before checking the 
level (see Fig 2) This delay is also frequency dependent and table- 
driven Note that if we simply delayed each cycle to keep from risking 



Dc inputs and low-level differential amplification for 
common-mode rejection are necessary features for an in- 
strument like the 8903A. One consequence of automatic 
ranging is that low-capacitance mechanical switching 
techniques cannot be used effectively. Needed are high- 
voltage reed relays, which affect high-frequency perfor- 
mance and require the use of compensation capacitors. 

For dc operation the first part of the input signal path is dc 
coupled with the input blocking capacitor bypassed. The 
output of the differential-to-single-ended amplifier can 
then bo monitored to obtain an accurately scaled represen- 
tation of the input dc level. 



Input voltages larger than 3V are attenuated by the input 
attenuator, a network of resistors that divide down the input 
signal. The appropriate tap point is selected by a reed relay 
network. If an overload occurs, the maximum attenuation 
setting is enabled. 

To protect the sensitive input amplifier following the 
attenuator from short-term transients, an overvoltage pro- 
tection network is used. For low-level signals the transfer 
impedance is low and signals applied to the input are 
coupled to the differential-input amplifier. However, for 
input signals large enough to damage the amplifier the 
output of the protection network is limited to a safe value. 



8 HFWLErr OAfKAHD JOURNAL AUGUST 19B0 



© Copr. 1949-1998 Hewlett-Packard Co. 



-Signal while settling 



Level looks good here 



Window of 
Acceptable Level 



V 



Time 

but slope measured in this interval 
shows need to delay to let signal settle 



Fig. 2. An example of slope checking by the 8903A micro- 
processor 

displaying invalid readings, every measurement would take longer, 
perhaps by several hundred milliseconds The slope checking gives 
us a rapid check on signal Quality which eliminates the need to delay 
every time 

After the slope check (and possible delay) we assume the output 
amplifier is well settled and ready for leveling. We then measure the 
voltage at its output, and if too high or low, ad|ust the gam accord- 
ingly After adjusting the gain, we make another pass through the 
leveling algorithm (unless we are now in the highest-gam range) To 
keep from getting caught in infinite loops, we require that after the first 
pass the gain never be reduced, and assume thai if il is, we have an 
unstable signal, which causes us to start the entire measurement 
cycle (tune source, level input, tune notch, level output) over again 
The final protection from infinite loops comes from counting the 
number of times we restart without displaying (each time we put a 
"— -" pattern on the display), and after 1 28 times we display Error 31 , 
which also goes out to the HP-IB 
Sweep 

Because the 8903A has the ability to generate accurate and com- 
paratively rapid measurements automatically, with no need lor the 
user to insert delays m HP- IB routines or wait for the display to settle 
before taking a reading, it is possible to have the instrument sweep 
itseil over a range of frequencies without user or computer interven- 
tion Here the calculating power of the microprocessor is brought to 
bear on the problem of determining the frequency increment lor each 
new point m the logarithmic sweep, based on the sweep range and 
the number of points per sweep (which can be set by the user). All the 
user need do to get a series of measurements spanning a frequency 
range is set the start and stop sweep frequencies and the number of 



points per sweep, and press the sago button This makes it easy to 
measure the frequency response of an amplifier 

X-Y Recorder Output 

The calculating power of the microprocessor-based controller .-s 
also apparent in the operation of the 8903A s X-Y recorder outputs 
Tnese outputs are driven by digitai-to-analog converters controlled 
by the microprocessor, rather than directly from internal detectors As 
a result, the user does not have to worry about the recorder output 
voltage abruptly changing when the analyzer autoranges The mi- 
croprocessor scales the recorder output according to the displayed 
reading and the plot limits entered via the keyboard. The recorder 
outputs are always between zero and ten volts so the recorder szero 
and vernier controls need be adjusted only once Thus, a properly 
scaled plot is easily generated by using the sweep and the X-Y 
outputs without any need tor an external controller. 

Special Functions 

Hidden behind the basic measurements are nearly forty special 
functions which provide extended measurement capability and 
many service aids For example, the analyzer can be given a load 
resistance in ohms and commanded to display ac level m watts 
Another special function changes the number of points per decade in 
a sweep, and several special functions modify display operation. 
Normally the left and right displays indicate the frequency and level 
(or distortion, etc ) of the signal applied to the analyzer input Some- 
times, when using the 8903A just as a source, the user may want the 
analyzer to display the frequency and level of the source Special 
function io provides this display 

Service aids provide front-panel display of many internal voltages 
and settings Without microprocessor control each special function 
would require one or more switches on the front panel instead of one 
special key. and would therefore probably not be included Thus, the 
processor allows implementation of useful features the user would 
not otherwise gel 
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The network consists of two back-to-back diodes which 
open up under large-signal conditions. 

The differential amplifier consists of three high- 
performance amplifiers. These have the necessary noise, dr. 
offset and frequency performance so that they do not de- 
grade the signal quality. The total input amplifier chain acts 
as a 4-dBstep amplifier with leveling hysteresis, which 
maintains the post-amplifier signal level within 6 dB (3 to 
1.5V rms). Should the output level change with time and 
deviate from this range, the gains of the attenuator, the 
diffen;ntial-to-single-ended amplifier, and the program- 
mable gain amplifier are adjusted to compensate. 



To summarize, the gain of the input amplifiers is mod- 
ified in three ways. First, the microprocessor monitors the 
input detector. If the detector voltage is too high or too low. 
the microprocessor varies the gain of the attenuator' 
amplifier chain to bring the level within bounds. Second, 
the overvoltage protection network limits if the input signal 
exceeds ±15V. Third, any sustained input overload trips 
the input overload detector. This detector monitors the 
input rms detector and the differential amplifier, and if 
either exceeds a certain absolute voltage limit, the overload 
detector trips, resetting the gain of the entire amplifier 
chain to its minimum value {maximum attenuation). This 
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Design for a Low-Distortion, Fast-Settling Source 



by George D. Pontis 



To fulfill the requirements of the 8903A Audio Analyzer the built-in 
source must have good performance in certain key areas First of all. 
tor swept measurements, it must be readily programmable and fast 
settling. It must also have low distortion and noise, and it must Mave 
very good amplitude accuracy over the entire frequency range of 20 
Hz to 100 kHz. 

The combination ol these conflicting requirements suggested the 
use of a high-performance RC oscillator instead ol a synthesizer The 
synthesizer is easy to program and settles quickly, but it is difficult to 
build a synthesizer with nose and distortion more than 80 dB oelow 
the fundamental Synthesized designs also do not have sufficient 
absolute level accuracy or flatness without leveling and a leveling 
loop that does not unduly degrade the distortion and settling lime 
would be very difficult to design 

Unfortunately, none of the common RC audio oscillator designs 
looked suitable either Usually the amount of distortion is inversely 
proportional to the settling time Also, the tuning elements usually 
float, making the circuit difficult to interface with programming lines. 

For these reasons a state-variable oscillator was chosen, similar to 
that proposed by Smith and Vannerson in 1975 ' Since this oscillator 
is built around a state-variable filter structure, inexpensive JFET 
switches can be used easily to switch the tuning elements More 
important is that the ALC design permits very rapid settling without 
trading off good distortion performance. Although the programmed 
frequency does not have the accuracy of a synthesizer, sufficient 
resolution is available to permit firmware tuning to within ±0.3% ot 
the programmed value 

This oscillator design can be described as a state-var.abie filter in 



which the resistor that determines the Q is replaced by an analog 
multiplier The control Signal for the multiplier is provided Oy an 
automatic level control (ALC) circuit The ALC circuit compares the 
oscillator amplitude with a stable dc reference obtained from a 
temperature-compensated reference diode The resulting error sig- 
nal is processed through two paths One path carries the cycle-by- 
cycle proportional error to the controlling multiplier This provides 
very last settling The other path includes an integrator in the loop to 
eliminate nearly all of tne steady-state error This design is theoreti- 
cally capable of settling the output amplitude within two cycles after 
small-signal amplitude disturbances 

There are two important refinements in the 8903A oscillator The 
first is the use of a special two-stage peak detector This consists of 
track,hold and sample/hold amplifiers to eliminate any aistortion- 
causing ripple on the detected peak output The second refinement is 
the addition of an ALC loop gam control to compensate the leveling 
loop, cycle by cycle, for changes in oscillator amplitude This greatly 
decreases the large-signal settling time of the oscillator, which is 
important when switching trom one range to another 

Fig 1 shows the oscillator integrators The gam constant of these 
integrators is changed in three-octave steps by selecting the feed- 
back capacitor. This gives us range switching. Coarse tuning within 
each range is done by switching the input resistors The eignt 
Dinary-weighted resistors provide 255 usable steps Placing the 
switching devices at the virtual gound point permits the use of JFETs 
with low drain-to-source on resistance (Ro S . un ). An individual transis- 
tor switch conducts when its gate is allowed to rise to ground poten- 
tial, and turns off when its gate is pulled to the negative supply, -1 5V 
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Fig. 1. The gam-switched inte- 
grator used in the 89034 's inter- 
nal oscillator Capacitors are 
switched to change ranges Resis- 
tor switching provides coarse tun- 
ing within each range 



reset occurs within ten milliseconds if the overload is se- 
vere. This protects the input from burning out or blowing a 
fuse, and allows for rapid overload recovery. 



Notch Filter 

The notch filter design challenge was twofold. First, it 
was necessary to design a low-distortion, low-noise filter 
that was also programmable. Second, to minimize overall 
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The switching scheme is economical because multiple resistor 
packages and quad comparator can be used to interface to TTL 
levels 

Fig 2 is a block diagram ot the oscillator integrators U1 and U2 
and inverter U4 form the state-variable filter structure Fine tuning is 
done by U4 This stage uses a switched resistor network similar to 
that used tor coarse tuning h the integrators As the resistors are 
switched the gam ot LM changes effectively altering the amount of 
signal transmitted from the output of U2 The gam is proportional to 
vA-Bfl, where R is me parallel combination of the selected input 
resistors and A and B are constants that provide a s5% fine tuning 
range This gives the oscillator enough resolution to tune within 
=0 2% of any frequency within the range ot the instrument 

It is generally true for sinusoidal oscillators that purity and settling 
time are chiefly functions of the ALC circuits or mechanisms used 
The two-stage peak-defector circuits are the key to the performance 
of the 8903A oscillator Oscillator amplitude data is obtained by the 
track'hoid and sample hold amplifiers in the following manner Switch 
S1 closes during the time the output is at its negative peak Capacitor 
C1 rapidly charges following the sme-wave amplitude up to its posi- 
tive peak At this time S1 opens, holding the peak voltage on C1 
Switch S2 then quickly closes and opens again, thus updating the 
sampled peak level held by capacitor C2. 

The two-stage scheme has several advantages For one. the first 
stage may be optimized for fasl daia acquisition, while the second 
state is optimized for long hold time Fasl acquisition is essential for 



good amplitude accuracy at high frequencies Low droop « impor- 
tant to maintain low distortion at tow frequencies Also this scheme 
has no steady-state ripple In the sampled output which would cause 
distortion to appear on the multiplier output and in turn on the mam 
oscillator output 

in practice, the two- stage scheme is also easie' to implement man 
a very fast smgie -stage sample hold circuit The required circuit 
functions are accomplished with simple JFtT switches and unity- 
gam buffers as shown in Fig 2 

For sampled data systems m general settling time is strongly 
dependent upon iocp gam For this Circuit the ideal integrator gam is 
linearly proportional to frequency Tnis effect is achieved by switch- 
ing the integrator resistor on once per cycle for a duration of 35 fis 
This process increases the integral error signal with frequency The 
duty cycle increases with irequency until the integrator switch S3 is 
left closed continuously for frequencies greater than 25 kHz Above 
25 kHz the oscillator easily settles in less than one millisecond 

8903A oscillator performance is largely limited by the quality of 
commercially available, reasonably priced analog multipliers The 
multiplier used is decoupled slightly (from optimum) to reduce its 
contribution to THD and noise This extends the osc illator settling time 
to a period of four to five cycles 

Reference 

1 E Varmeison anil K Jmn'i A Lc* D-alamo" OiCiHaw With 'a r ,\ MApHtfft SlaObja- 
lion lnie'i-OK>nal Jou'nw si Elejlro'xs W 39 No 4 '975 oo 465-172 




ALC Output 



Fig. 2. Block diagram ol Ihe 8903A s internal stale-variable oscillator 



measurement time, it was necessary to develop an accurate, 
rapid fine-tune mechanism with quick recovery from over- 
loads and mistiming. To achieve good (J over the frequency 
range, an active KC filter is necessary. To tune this device, a 



variable resistance or conductance device is needed, since 
achieving the tuning range with variable capacitors or in- 
ductors is impractical. Many resistively tunable active con- 
figurations arc feasible. However, determining the op- 
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Floating a Source Output 



by George D. Pontis 



To provide the greatest versatility in both benchtop and systems 
applications, the 8903A Audio Analyzer s built-m source is tloating. 
This lets the user eliminate ground-loop errors, sum signals, and add 
dc offsets to the source output. 

Previous designs used a separate, isolated power supply for the 
source circuits- This method is straightforward and offers very good 
low-frequency common-mode reiecnon. However, there are several 
reasons why this arrangement is not used in the 8903A. 

The biggest problem with the floating power supply approach is 
interlacing with the digital programming lines Since only three of the 
thirteen attenuator lines have relay isolation and none of the nineteen 
oscillator lines can be floated, over thirty lines must be coupled in 
some manner. One solution is to float only the final output stage. This 
eliminates the need tor couplers, but requires a high-performance 
differential-input amplifier to reiect the common-mode signal that 
appears at the input of the floating stage. Since a floating power 
supply is still required, the cost of this approach is relatively high 

The 8903A solves this problem with a single-ended-to-differential 
output converter This circuit, shown in Fig. 1 , operates on the instru- 
ment s ground-referenced ± 15V supplies and requires only two op- 
erational amplifiers A precise combination of negative feedback, 
positive feedback, and cross-coupling yields a symmetrical differen- 
tial output with infinite common-mode reieclion and a well-defined 
output impedance 

An analysis of this circuil is generally a tedious procedure because 
of the number of components involved. However, me high degree of 
symmetry in the circuit can be exploited to great advantage by using 
the relations R2/R1 = R12/R11, R3 = R7 R6 = R10, and R4 = R5 = 
R8 = R9. From these relationships, one can derive the expression 
R2'R1 = (2R6 + R4).'2R3, which is a necessary condition for achiev- 
ing an infinite common-mode output impedance. Then it is easy to 
calculate the differential output impedance and the open-circuit volt- 
age gam. The resulting equations can be manipulated to find suitable 
values. For the resistor values used in the 8903A, the associated 
gam is 1 125 and the output impedance is 480 ohms. The output 
is further padded with a 1 20-ohm resistor to yield the desired 600-ohm 
output impedance. 

It would have been possible to use resistors that gave an output 
impedance of exactly 600 ohms instead of 480 ohms but this would 
have required setting up and stocking a supply of several extra odd 
resistor values. As it is, the circuit is realized using 0.1%. 25-ppm 
resistors that are also used elsewhere In the instrument 
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Fig. 1 . Smgle-ended-lo-dillerential output converter provides 
a lloatmg output lor the 8903 As interna/ oscillator 

The easiest way to see how this circuit works is to eliminate either 
the inverting (low) or noninverting (high) half of the circuit by shorting 
the respective output to ground. Fig. 2 shows the reduced circuit 
when the low half is grounded. If R4 is disconnected, the circuit will 
have a forward gain of about one, and an output impedance of 276 
ohms. R4 works in conjunction with R6 to provide voltage and current 
feedback that causes the gain and the output impedance to rise. 

To demonstrate that the output is truly floating, we ground the input 
and apply a test source to both outputs. Ideally, the current flow from 
the test source should be zero. Fig 3 shows a block diagram and the 
reduced circuit for this test Here it can be quickly calculated that the 
output of U1 will rise just enough over that of the test source to make 
the current through R6 cancel the current through R4 and R8. Note 
that the current flowing through sources V1 and V2 is supplied by the 
other half of the circuit, which is not shown 



timum match between the filter configuration and the vari- 
able resistive element is not straightforward. 

Let's go through the alternatives and the tradeoffs, 
Switchable resistor networks have good distortion and 
noise characteristics but do not provide continuous tuning 
coverage and require extensive switching circuitry. Photo- 
resistors can be driven over a large resistance range and 
provide continous tuning. However, the noise and distor- 
tion they add to a signal are greater than the required level 
of 90 dB below the signal level. They can be used as fine- 
tuning elements if coupled only partially into the circuit. 
These devices can also be slow and are awkward to control 
rapidly, reducing the tuning speed. Finally, they tend to 
vary significantly from device to device and with time and 
temperature, making compensation difficult. 



Four-quadrant analog multipliers also do not have the 
90-dB performance necessary, but they too can be lightly 
coupled into the circuit for fine tuning. These devices are 
fast, inexpensive, and easy to drive. There are many varia- 
tions on this type of circuit, some of which can be obtained 
in integrated form. Those most suitable use a differential 
pair of bipolar transistors as a variable gain element by 
varying the common-mode current. 

Light bulbs as variable resistive devices are relatively 
linear but are slow and have a limited dynamic, range, 
Thermistors, diodes and other nonlinear devices would all 
be useful only for fine-tune applications. The drive and 
compensation circuitry for all of these alternatives would 
be complex and the overall performance marginal. 

The tuning elements selected were switchable resistor 
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Fig. 2. Eliminating the inverting side of the circuit ot Fig 1 by 
shorting the low output to ground results in this reduced cir- 
cuit 

In practice it was found that parasitic effects and the 
electromagnetic-interference (EMI) filters degrade the circuit bal- 
ance when the frequency approaches 100 kHz However, each 
board is tested for a minimum of 50 dB common-mode reaction at 1 
kHz A typical unit has greater than 40 dBreiection at 100 kHz Also, a 
1 0-kfl resistor internally ties the low output to the chassis ground. This 
provides a reference for the output when no external load is con- 
nected. Without this resistor, the common-mode output voltage is 
indeterminate. 

One initial concern about the circuit was difficulty of troubleshoot- 
ing problems, such as one of the resistors drifting out of tolerance, 
causing poor common-mode reaction This problem was solved by 
implementing the following test procedure. First, the input to the 
floating amplifier is set to exactly 1 00V rms using the special func- 
tions built into the 8903A. Then, the technician shons the low side to 
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(b) 

Fig. 3. la} To demonstrate that the output of We circuit is 
floating, the input is grounded and a voltage source is applied 
to both outputs . (b ) Reduced circuit tor this test (only hall of the 
circuit is shown). Current through R6 cancels the current 
through RA and R8 Current through VI and V2 is supplied by 
the other hall of the circuit (not shown) Thus no current is 
drawn from the test source 

ground and measures the potentials al several circuit nodes The 
measured values can be compared to the calculated values pub- 
lished in a table in the 8903A manual A second set of measurements 
can be made, if necessary, with the high side grounded By observ- 
ing the deviations between the measured and calculated values, it is 
easy to locate the faulty component 



networks for low-distortion coarse tuning, and a four- 
quadrant multiplier for fine tuning between the discrete 
steps of the switchable resistors. Since the four-quadrant 
multiplier is coupled into the circuit only enough for ±7% 
tuning, it does not contribute significantly to the overall 
noise and distortion. A resistor switching network is best 
implemented if one end of the network is dynamically kept 
at ground potential: this relieves many constraints on the 
switching network. To this end. a state-variable* notch 
configuration is used (see page 14). With this design, low- 
distortion tuning over a 10:1 frequency range is achieved. 

"The namo 'state-ya'iable titte'' >s 3 consequence of the laci that Ihe equations mat » 
scribe the system can be wrilten in a lorm thai MS Ihe classical stain-space description ol 
a lineal system, i e * - Ax I Bu What* » is the response lor stale variable) ol the system 
10 an input u 



Capacitors are switched into the network to change fre- 
quency in three-octave bands and provide complete cover- 
age of the frequency range of the analyzer. 

To complete the fine-tuning path, a synchronous detector 
mixes the filter output waveform with the fundamental 
waveform. If any fundamental exists in the notch output, a 
dc current is generated to fine tune the notch. The critical 
parameters here are 100-dH dynamic rangeand rapid opera- 
tion. A PET double-balanced mixer was selected. The input 
mixing signal that drives the FET is a square wave, rich in 
odd harmonics, so the circuit responds to odd-order har- 
monics as well as to the fundamental. This is the classical 
solution. A complete null may not hi; achieved if third, fifth, 
or higher-order odd harmonics are present. Total error can 
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A Digitally Tuned Notch Filter 



by Chung Y. Lau 



The notch filler in ihe 8903A Audio Analyzer rejects the fundamen- 
tal frequency component of the incoming signal. This filter consists of 
a state-variable active filter and fine-tune and fine-balance control 
circuits. A simplified schematic of the filter minus the control circuitry 
is shown in Fig 1 . 

Conventional notch filters (bndged-T, etc ) used in distortion 
analyzers are not well suited for digital control because they require 
expensive relays or analog switches The state-variable filter ap- 
proach has the following advantages 
1 The Q of the filter is fixed and independent of frequency 

■ Tuning is accomplished by switching resistors and capacitors in 
the two integrators. Inexpensive JFETs are used as switching 
elements and the switch drivers are simple because both the gate 
and the channel of the JFET are virtually at ground potential when 
the switch is on. 

■ Three filter outputs (low-pass, bandpass, and high-pass) are 
available simultaneously The low-pass and high-pass outputs are 
phase-shifted from the input by 90° when the filter is tuned to the 
input frequency. These two outputs are used in the control circuits, 
so extra phase-shitting circuits are not needed 

- Distortion generated in the fine-tune circuitry is filtered by two 
integrators before appearing at ihe notch output 
In distortion mode, the filter is tuned in the following manner. The 
microprocessor counts the input frequency and tunes the notch filter 
to the same frequency by switching in tne proper capacitors and 



resistors When the filter is tuned to the fundamental frequency of the 
input signal the bandpass output (\4 P in Fig. 1) inverts the funda- 
mental component of the input signal U5 sums V| N and V BP to cancel 
the fundamental component of V IN exactly However, because of the 
phase and amplitude characteristics of V BP , the harmonics present in 
V, N are relatively unattenuated. 
Analytically, we have 



V BP 
— (s) = 



-sw,,'Q 



s' + stu'Q + ai- 



V BP + V lt 



•<s) = 



S 2 + ^ 



"IN 



S z + S<u 0 'Q + ml 



where oi Cl is the center frequency in rad.'s. Therefore, there is zero 
transmission at s = ±\m 0 . 

The state-variable filter alone can provide only about 15 dB of 
fundamental reaction because of the coarseness of the frequency 
tuning (use of discrete resistor values, resistor mismatches, etc.). 
Fine-tune and fine-balance control circuits are necessary to achieve 
a notch depth greater than 90 dB. The fine-tune circuit insures that the 
filter is exactly tuned to the input frequency (no phase error) and the 



Correction Signals 
from Control Circuits 




Fig. 1 . 8903A notch filter. Correc- 
tion signals from the control cir- 
cuits ol Fig. 2 provide line tuning 
and tine balance 



be as much as 0.46 dB, but this is considered reasonable. 

Rapid response is achieved in two ways. First, when an 
overload occurs, the circuit gain is quickly reduced by the 
output amplifier circuit, minimizing the transient and thus 
minimizing the recovery time of the synchronous detector 
from the overload. Second, the response time of each fre- 
quency band has been optimized. Theoretically, it is impos- 
sible to fine tune a filter as rapidly at 20 Hz as at 1 kHz. It 
simply takes a much longer time to detect a null at 20 Hz, 
but even more important is that the time response of a 20-Hz 
notch filter is much greater. In fact, a fundamental burst 



tone applied at the input of the notch will not immediately 
be nulled at the output. The entire input will first appear at 
the output and then decay away in proportion to f(notch) - Q 
where Q is the Q of the notch. Thus, as the notch frequency 
is increased, so is the speed of the fine-tuning circuitry. The 
bottom frequency band (20-200 Hz) is relatively slow but 
nulls optimally at 20 Hz. For higher frequencies the instru- 
ment response time improves. 

Output Amplifier 

For the amplifier following the notch filter, gain accuracy 



14 HEWLETT-PACKARD JOURNAL AUGUST 1980 



© Copr. 1949-1998 Hewlett-Packard Co. 



*15V 




Fig. 2. 8903-4 notch filler line tuning and balance circuits. 



fine-balance circuit insures thai |V a p/V| N | = 1 at the input frequency 
(no amplitude error). The fine-tune circuit is described in some detail 
here; the fine-balance circuit operates in a similar manner. Fig 2 is a 
simplified schematic of both circuits. 

The fine-tune circuit operates as follows (see Fig 2). The low-pass 
filter voltage V lP , which is phase-shifted 90° from the input, drives the 
tune comparator U6A, which turns JFET switch Q5 on and off When 
Q5 is on. point A is essentially grounded and no current Hows into the 
integrator capacitor C3 When Q5 is off. the notch amplifier feeds the 
tune integrator Because of this chopping action, the dc current that 
(lows into the integrator is caused by the notch output component 
that is synchronous with V LP .* 

The tune integrator output is a dc voltage that can be changed only 
by the dc current flowing into the integrator. This voltage feeds one 
input of the multiplier U9A The other input to the multiplier is V BP The 
product of the two inputs is summed into the state-variable filter via 
U3 (see Fig 1). The net result is that the tune integrator voltage can 
change the effective resistance (R t in Fig 1) and hence the center 
frequency of the notch filter 

The direction of change is such that any notch output components 
in synchronism with V lP are reduced At steady state, the 

■Actually third harfnomcs <n Ine notch output can aiso cause dc current to How in trie 
integrator This error causes a small amount ot the lunaamenta! frequency component to 
pass through the notch litter However tne maximum error contribution to distortion mea- 
surements Irom this ettecl is only aoout 0 46 dB 



tune voltage is stable, although not necessarily zero, and no dc 
current flows into the tune integrator Therefore, there is no notch 
output component synchronous with V LP . 

The fine-balance circuit insures that there are no components 
synchronous with V BP in the notch output V BP and V LP are both at the 
fundamental input frequency and are in quadrature with each other 
Thus, at steady state there is no fundamental frequency in the notch 
output 
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and good frequency response are important. In ac level 
measurements the signal travels through the output 
amplifier chain and is detected by the output rms detector. 
Any error degrades the 8903A's performance. The output 
amplifier allows the low-level distortion products leaving 
the notch to be accurately detected by limiting the required 
dynamic range of the rtns-to-dc converter. The rms detector 
is accurate over only a 30-dB range, and the output 
amplifier boosts these signals into the range of the detector. 
As with the input circuitry, rapid recovery from overload 
conditions is crucial. If the notch becomes mistuned be- 
cause of a disturbance at the input, the output suddenly 



increases dramatically, sending the amplifier into overload. 
This in turn generates dc offsets in the amplifier chain that 
can take seconds to decay, even if the amplifier may start 
operating again much sooner. Thus a large low- frequency 
impulse appears at the amplifier output along with the 
signal being amplified. The composite signal is transferred 
to the output rms detector, which responds to a time- 
weighted average of the total. II the impulse is a significant 
fraction of the signal, the rms detector will not give a true 
indication of the signal amplitude. This can cause problems 
when the instrument ranges automatically. In effect, it 
forces the leveling algorithm to wait much longer to con- 
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firm that the amplifier output is within the detector's range. 
This in turn slows down the rate at which the instrument 
can determine the proper measurement range and display a 
reading. Even if the impulse is 20 dB less than the signal, 
the detector error can be as much as 0.5%. This will also 
increase the amount of time required by the instrument to 
make an accurate reading from the detector once the proper 
range has been obtained. 

To alleviate these problems, the size of the transient is 
minimized in three ways. First, operational amplifiers and 
circuit configurations are used that have better than average 
overload immunity. Second, an overload detector is placed 
at the output of the rms detector. If the signal level becomes 
too large, the overload detector trips and the amplifier gain 
is reset to 0 dB. Third, a 1 3-1 1 z high-pass filter is placed 
before the output detector. This significantly reduces the 
duration and amplitude of any transient and hence keeps 
the transient from significantly increasing the total mea- 
surement time. The only delay factors that remain are the 
controllable and predictable settling times of the notch 
circuit and the rms detector. 

The response time of the output rms detector is a com- 
promise between rapid settling and low-frequency accu- 
racy. A configuration was selected that settles to within 1% 
of a 10:1 step in 350 milliseconds, and has a steady-state 
error of 0.2% at 20 Hz. This includes filtering in the detector 
and additional filtering following the detector to reduce 
excess ripple. For leveling purposes the ripple is not sig- 
nificant, so the microprocessor uses the detector output 
when leveling and avoids the extra delay contributed by the 
additional filter. The output detector and filters could have 
been designed with switchable time constants to respond 
more rapidly for higher-frequency signals. However, the 
penalties would have been additional circuit complexity 
and the ambiguity of not knowing when to invoke the longer 
time constant. A 20-kHz signal, for example, might still 
have a significant low-frequency component, which would 
cause excessive error with a more rapid time constant. 

Oscillator 

Many of the design considerations for the notch filter 
apply equally to the oscillator. In both cases tuning consid- 
erations were the same, with switchable resistor networks 
used as the decade tuning elements and four-quadrant 
multipliers used for amplitude control. In many ways the 
oscillator and notch circuits can be seen as duals. The notch 
generates a pair of zeros on the jtu axis that reject the funda- 
mental component, while the oscillator generates a pair of 
poles on the jto axis that generate sustained oscillations. The 
trick in the oscillator is to keep the poles exactly on the axis 
to maintain constant output amplitude. This must be done 
continuously by the automatic leveling circuit. If the fre- 
quency of the circuit deviates from the desired frequency, 
the circuit can be fine tuned by the microprocessor, which 
monitors the output frequency on a sampled basis. The 
major performance goals of the oscillator were low noise 
and distortion, rapid amplitude and frequency settling, and 
digitally programmable frequency control. Again the 
state-variable filter configuration along with a special level- 
ing circuit offered the flexibility and performance required. 
The oscillator design is described on page 10. 
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It was determined during development that the oscillator 
would have to run at a constant output level to maintain 
reasonable settling and noise performance. It was also de- 
sired to have a floating output. The attenuator and output 
amplifier circuit (see page 12) takes the oscillator output 
level and translates it to the selected floating output 
amplitude. To minimize cost and achieve overall output 
accuracy goals the attenuation is done in two stages. Coarse 
amplitude steps are implemented with a 2.5-dB/step at- 
tenuator network. Smaller steps are provided by a resistive 
ladder network that adjusts the amplitude linearly in small 
discrete steps. The combination can adjust the amplitude 
within a nominal *0. 15% worst case. Computation of the 
proper switch settings is an easy job for the computational 
skills of the microprocessor, 

EMI Design 

Meeting the required electromagnetic interference [EMI] 
and susceptibility goals was a bit more challenging than 
initially expected. Large-amplitude RF fields tend to gener- 
ate voltages on exposed cabl ing and circuits. These voltages 
overdrive many of the active circuits, causing nonlinear 
operation and distortion. To avoid direct exposure to these 
fields, the analog circuits are housed in an internal EMI- 
tight box. The box has an aluminum frame around the sides. 
The bottom cover is the ground plane of a printed circuit 
board and the top cover is a removable EMI-tight lid. Re- 
moval of the lid. which is held in place by only two screws, 
makes all the circuits available for service. The micro- 
processor boards are sufficiently shielded by the instru- 
ment cabinet and do not require the extra shielding. To keep 
the RF fields from developing voltages on the cabling feed- 
ing the circuits, special precautions were taken. First, from 
the inner box to the front panel, shielded cable is used. 
Second, BNC connectors are provided on the front panel. 
The BNC connectors allow the attachment of shielded ca- 
bles directly to the instrument if desired, thus preventing 
EMI pickup. The instrument's digital circuitry also gener- 
ates EMI related to harmonics of its 2-MHz clock. This 
problem was minimized by means of RF gaskets on some of 
the cabinet seams and by installing an EMI suppressing 
filter on the power line input. As a result, the instrument 
will not disturb sensitive receivers operating nearby, and 
yet will perform well near a powerful transmitter. 

Frequency Measurement 

A key feature of the 8903A is its ability to measure fre- 
quency automatically, even when the input waveform may 
have a significant amount of noise and distortion and the 
amplitude may vary from 6 mV to 300V. Part of this problem 
is solved because the instrument is autoranging and keeps 
the leveled waveform within 6 dB over most of the input 
amplitude range. But before the signal can be accurately 
counted it must first be converted into a binary signal hav- 
ing the same period as the major frequency component in 
the waveform, and herein lies a problem. If a zero-crossing 
circuit is used, noise may cause multiple crossings and a 
false indication of the frequency. Hysteresis in the detector 
will help, but if the hysteresis is too large, smaller- 
amplitude waveforms may not trigger the detector at all 
while large-amplitude waveforms will have relatively little 
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hysteresis protection when large noise components are 
present- To alleviate this problem, the 8903A employs vari- 
able hysteresis. As the peak amplitude of the signal varies, 
so does the hysteresis level, which is maintained at approx- 
imately one- half the positive peak for the positive portion of 
the waveform and approximately one-half the negative 
peak for the negative portion of the waveform. Hence noise 
immunity remains constant regardless of the incoming 
waveform. Hysteresis is implemented with a bipolar peak 
detector followed by a dual comparator. The waveform is 
transferred to a reciprocal counter, which measures the 
period of the signal, and the microprocessor inverts this 
period to get frequency. 
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A Custom LSI Approach to a Personal 
Computer 

by Todd R. Lynch 



THE NEW HP-85 PERSONAL COMPUTER, which 
was featured in last month's issue, is a system totally 
integrated into a single package. Included in this 
system are a CRT. printer, tape drive, and keyboard. To 
control the I/O (input/output) devices and to interface to 
ROM (read-only memory) and RAM (random-access mem- 
ory), the design uses LSI (large-scale integrated) circuits 
designed and fabricated by Hewlett-Packard. 

The design of custom chips minimizes the cost of the 
electronics. Also, the power dissipation is reduced to a level 
that permits the use of air-convection cooling, eliminating 
the need for a fan. The LSI designs save large amounts of 
printed-circuit-board space, making a small system pack- 
age possible. By designing LSI circuits dedicated to each 
mechanical subassembly, features can be added to the over- 
all system that would be nearly impossible to incorporate 
with discrete logic designs. 

HP-85 LSI Chip Set 

There are nine custom LSI circuits in the HP-85. These 
circuits are interconnected as shown in Fig. 1. Eight of the 
LSI circuit designs use NMOS technology. The ninth circuit 
uses a bipolar technology with two layers of metallization. 

The heavy emphasis on LSI design makes the HP-85 very 
compact. The electronics in the machine consists of these 
circuits plus a power supply, clock generator. CRT raster 
scan circuitry, and various motor and printhead drive cir- 
cuits. 

The system is partitioned as follows. 
- CPU (central processing unit). The CPU commands the 
bus control lines and interfaces directly to the system 



operating commands (firmware) stored in the four ROMs. 
ROMs. Each ROM chip has the same basic custom design, 
differing only in the bit pattern permanently stored on 
each chip. 

■ Read/Write Memory. The user memory consists of eight 
commercially-available dynamic 16Kxi RAMs. 
RAM Controller. This chip is designed to interface be- 
tween the CPU and the RAMs. 

Buffer. This chip is designed to provide the capability for 
system expansion by adding plug-in units in the rear of 
the machine. 
Keyboard Controller. 
Printer Controller. 
CRT Controller. 

Display Memory. Four commercial dynamic 16Kxi 
RAMs are used to store data for the display. Approxi- 
mately one-fourth of this memory is used for the alpha- 
numeric data and the remainder is used for graphic data. 
Cartridge Controller. 

Sense Amplifier. This bipolar chip interfaces the tape 

cartridge controller to the magnetic tape head. 

Table I is a summary of the custom LSI chip sizes, number 
of pins (wiring connection points), and power dissipation. 

The machine's sixteen-bit address allows direct access to 
65.53(5 bytes of information. The memory map in Fig. 2 
shows how this space is allocated. Note that 32K bytes of the 
address space are devoted to ROM while lfiK bytes are for 
RAM. Additional RAM and ROM capacity may be added to 
the system. The I/O address space occupies the upper 256 
bytes of memory. Each I/O device controller has from two to 
four dedicated addresses assigned to it. 
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Fig. 1. The HP-85 system block 
diagram contains nine custom LSI 
circuit designs The user and dis- 
play memories are the only parts in 
the diagram that use commercially 
available circuits. The develop- 
ment of custom circuits enabled 
the system to be contained in a 
compact package at low cost and 
without the need tor a cooling tan. 
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Table I 

Custom LSI circuit characteristics All of the circuits are fabricated 
with NMOS silicon gate technology except for the sense amplifier 
which is made with dual-layer-metallization bipolar technology 



CPU Design 

Several commercially available CPUs were considered at 
the beginning of the project, but none could provide all the 
features needed to efficiently implement a powerful scien- 
tific BASIC language machine. BASIC requires high- 
precision arithmetic, which is best accomplished with dec- 
imal rather than binary numbers. Many different stacks are 
needed to parse (separate a statement into executable steps) 
and execute the language. The ability to handle variable- 
length data is required for variable-length tokens (bit se- 
quences from one to several bytes in length), and multilevel 
vectored interrupts are needed to handle the I/O devices in 
real time. The design of the custom NMOS CPU incorpo- 
rates all of these requirements plus many more. 

Fig. 3 shows a block diagram of the HP-85 CPU, Major 
blocks are the 7000-bit PLA (programmable logic array), 
64-byte register bank, and eight-bit ALL' (arithmetic logic 
unit) and shifter. Each CPU instruction is decoded by a 
microprogram in the PLA that directs the rest of the chip to 
perform the desired function. 

A very powerful feature of the CPU, which is incorpo- 
rated into the Pl.A microprogram, is the ability to handle 
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Fig. 2. The HP-85 memory map allocates lour 8K byte system 
ROMs and one 16K-byte internal RAM The RAM is used lor 
user memory and 256 dedicated HO locations By using bank 
selection, up to six more 8K-byte ROMs may be added. The 
RAM may be expanded by 16K bytes it a plug-in module is 
added to the system. 
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Fig. 3. Trie CPU in the HP-85 is a custom NMOS circuit de- 
sign A 7K-bit programmable logic array (PLA) controls the 
interactions between the register bank, the ALU. and the rest 
ol the HP-85 system. 

data ranging from one to eight bytes in length. This multi- 
byte feature lets the programmer, for example, add two 
eight-bvte mantissas, increment a sixteen-bit address, or 
load into the CPU a three-byte token from memory, each 
with a single instruction. With an off-the-shelf CPU. this 
would require setup and iteration within a software loop. 

The eight-bit ALU is capable of shifting and can do both 
decimal and binary arithmetic. The programmer sets a 
mode bit to specify the type of shifting or arithmetic de- 
sired. This, combined with the multibyte feature, lets the 
programmer easily work with signed, floating-point man- 
tissas up to sixteen digits in length or two's-complement 
binary integers up to sixty-four bits in length. 

Of the 64 eight-bit registers contained in the CPU. one 
pair is dedicated to the program counter, one pair to the 
stack pointer, and one pair to internal index calculations. 
The rest are general-purpose registers. One advantage of 
having a register as the program counter is that it can easily 
serve as one of the operands for any CPU instruction. If it is 
modified by an instruction, then an immediate jump to the 
new location occurs. The stack pointer points to the sub- 
routine return address stack stored in memory. Additional 
stacks can be created with any other consecutive pair of 
registers in the CPU. Thus, the programmer can maintain 
many different stacks at any one time. Sixteen instructions 
are dedicated to manipulating data on the currently desig- 
nated stack. Any pair of consecutive registers can also be 
used as a sixteen-bit index register. This lets the program- 
mer index into several data arrays at the same time. 

To handle real-time I/O devices, the CPU has multilevel 
vectored interrupt handling ability. Up to 127 interrupt 
vectors can be accommodated by the architecture. The CPU 
can also be halted by an external device. This lets that 
device control the system bus. so it can have direct memory 
access at high speeds. 

No accumulator is present in the CPU. This is made 
possible by the design of the (>4-byte register bank. The 
registers are constructed as a two-read/one-write memory. 
This means that, at a given time, any two bytes can be read 
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Q Bit Line 



O Bit Line 




Fig. 4. Independent control ol each word line is used to 
achieve a cell that can be readout onto either bit line. The RAM 
made up with these cells can be read two words at a time. 



from the memory and operated on in the ALU. The result is 
returned to one of the accessed byte locations. This se- 
quence takes one processor cycle (1.6 £is). 

To design a two-read/one- write memory, a standard static 
RAM cell with a special variation is used (see Fig. 4). In a 
normal static RAM cell, a common word line enables both 
transmission gates between the bit lines and the cell. In a 
two-read cell, there must be two different word lines. Word 
line A enables complement information from one cell to be 
read out onto the Q-bit line while word line B enables true 
information from another or the same cell to be read out on 
the Q-bit line. 

The circuit design of such a cell must be done carefully. 
Since the bit lines are precharged before reading, it is pos- 
sible that by enabling a single word line, a cell could be 
made to flip (change state! rather than be read. Inadvertent 
flipping is prevented by using a low voltage for the word 
line and a proper size ratio between the cell's pull-down 
transistor and the transmission gate. 

The static RAM cell's pull-up device also received close 
scrutiny during the design phase. A cell with a low-power 
depletion-load pull-up requires an area larger than 6400 
square micrometres (10 square mils) and a quiescent power 
of 120 microwatts. By using a polysilicon pull-up resistor, 
the RAM cell area is reduced by 40% and the power is 
reduced by 80%. The savings in area amounts to a reduction 
in size of more than 0.25 mm on each side of the die. The 
reduction in total power consumption is fifty milliwatts. 

To obtain these reductions, a polysilicon resistor process 
was developed for the CPU. An additional masking step is 
required to define the regions on the chip where the 
polysilicon layer is lightly doped, hence creating the poly- 
silicon resistors. The doping level was chosen for a sheet 
resistance of 10 7 ft/a. An undoped layer could not be used 
because the sheet resistance was so high that the junction 
leakage currents at elevated temperatures resulted in un- 
welcome voltage drops across some resistors. It was also 
discovered that if the contact mask overlapped the polysili- 
con resistor area, the aluminum metallization could spread 
into the lightly doped resistor and lower its resistance. 



System Control and Timing 

Eight bus and three control lines leave the CPU. The 
eight-bit bus is used to time-multiplex a sixteen-bit address, 
instructions, and multibyte data quantities. The three con- 
trol lines — LMA (load memory address), RD (read) and WR 
(write) — indicate to the system circuits what type of infor- 
mation is on the bus. 

When LMA is low, all chips in the system know that one 
byte of a two-byte address will be placed on the bus. The 
LMA signals always come in pairs since all addresses are 
sixteen bits. Each chip in the system reads this address and 
compares it to the address range to which the chip is sup- 
posed to respond. Setting RD low indicates that the CPU 
wants to read the contents of the address most recently sent 
out. When WR goes low. the CPU wants to write into the 
location most recently addressed. 

The multibyte feature is accomplished by sending out an 
address with two LMA signals followed by one to eight RDor 
WR signals. The circuit being read from, or written to, is 
expected to increment its memory address register every 
time it sees a RD or WR. The CPU sometimes can fetch 
consecutive instructions from memory by simply sending 
out additional RD commands. This speeds up the instruc- 
tion rate of the machine since sending out an address for 
every instruction and piece of data is not necessary. 

Four nonoverlapping phases (Pig. 5) with 200 ns width 
and 200 ns spacing clock the circuits in the HP-85. In one 
cycle the system can access a preaddressed memory loca- 
tion, read it into the CPU, add it to a CPU register and store 
the result in the CPU register. An eight-byte add requires 
eight cycles plus the time to fetch the command, which is 
usually three cycles. 

Before any chip writes on the bus, the bus gets precharged 
to a logic one. Therefore, the circuit desiring to put a one on 
the bus must merely continue to hold the bus high. If it 
wants a zero, it must discharge the bus. In NMOS technol- 
ogy it is easier to discharge a bus than to bring it high. 
Consequently, the circuits were designed with large de- 
vices to discharge the bus and relatively small devices to 
maintain a high level on the bus. This minimizes the chip 
area required for each circuit's driver logic. 

The svstem bus has fewer spurious transitions using the 

TL 



<412 

..1.2 

rf,21 
LMA. 55. wfi 

Bus to CPU 



Valid 



Precharge 



Float Precharge 



Valid 



Bus from CPU Float 



Precharge 



Valid 



Fig. 5. HP-85 system timing and control The control lines are 
valid by the time the <4' 2 clock pulse occurs. II a system chip is 
to be read, it must respond during the tb2 clock pulse. Informa- 
tion then enters the CPU, goes through the ALU, and is stored 
during the next $1 clock pulse. When the CPU sends out 
addresses or data, they are valid by the time the <t>1 clock 
pulse occurs. 
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precharged scheme. Timing is such that the driving circuit 
is not enabled until its data is valid and ready to go on the 
bus. It is not possible to discharge the bus inadvertently 
when the chip wants a logic one to be output. 

ROM 

The NMOS system ROM is organized as an 8K by eight- 
bit array. Because this circuit must respond to multibyte 
transfers, its memory address register is designed to incre- 
ment as RD commands are given. The circuit can be enabled 
or disabled by a bank select command. When plug-in ROMs 
are used, it becomes necessary to selectively enable or dis- 
able ROMs in the address space from 24K to 32K (Fig. 2). 
Each ROM in this address space has a unique eight-bit bank 
number, and only one bank can be active at any given time. 

The core of the ROM consists of 64-input. minimum 
geometry. NAND gate arrays. The principle of operation is 
to precharge the NAND stack from both ends. This requires 
600 ns. At the beginning of phase <b12. one 64-device stack 
per bit of output is selected. If the stack discharges, the bus 
bit will be a zero, otherwise it will be a one. The transistors 



in the stack are preprogrammed during fabrication with the 
depletion mask. A depletion transistor is a logic zero, an 
enhancement transistor is a logic one. A high voltage on 
either transistor causes conduction, while a low voltage 
causes conduction only through the depletion device. 

When trying to discharge a 64-device stack. 63 inputs are 
high and the selected bit is low. If a depletion device is at 
that location, the stack discharges, producing a zero. 

RAM Controller 

The NMOS RAM controller chip interfaces the CPU to 
eight 1 6K dynamic RAMs. Because it is a memory controller 
chip it has an incrementing memory address register. Also, 
the address space to which it responds is mask program- 
mable. Therefore, a RAM controller chip different from the 
one in the mainframe of the HP-85 is used for the 16K RAM 
plug-in. This chip knows not to respond to dedicated LO 
addresses in the upper 256 bytes of memory. 

An important function of the RAM controller chip is 
refreshing the dynamic memory. An internal timer on the 
chip tells it when the next sequential location is due to be 



The HP-85 Software Development System 

by Nelson A. Mills 



The HP-85 is based upon a new eight-bit custom processor One of 
the first steps in the development of the machine was to provide a set 
of software development tools. The HP-85 software development 
system consists of an assembler, a hardware interface for the HP-85 
simulator, and software debug system (see Fig. 1) 

The assembler for the HP-85 processor was designed to run on an 
HP 1000 Computer The assembler supports the full range of the 
HP-85 processor's instruction set and provides several useful 
pseudo-operations Programs may be either absolute or relocatable 
and the program origin may be reset at any time Several instructions 
are provided to facilitate data definition and may exist locally within 
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HP-85 Eleclronics 

Fig. 1. HP-85 software development system 



the program, or may be retrieved from a file of global data definitions. 

The hardware interface designed for the HP-85 software develop- 
ment system provides the ability to use an HP 1000 to control execu- 
tion of the HP-85 simulator The interface contains an on-board RAM 
which can be downloaded from the HP 1000 and is used to simulate 
the HP-85 ROM Included are two breakpoint registers which can be 
used to halt execution of the HP-85 processor at either of two 
specified addresses. The processor can be made to execute steps in 
either a continuous-run mode, or in a single-step mode where execu- 
tion halts after each instruction The interface also provides the ability 
to halt the HP-85 processor at any time, and to read data from or write 
data to the system RAM or CPU registers 

The software debug system, provided as part of the development 
system, runs in the HP 1 000 and performs four important funciions. It 
includes a relocatable loader that is used to download the RAM on the 
interface board with the software under development It controls 
execution of the HP-85 simulator via the interface board. It dispiays 
the current status of the breadboard, including CPU status, program 
counter, and the current contents of all CPU registers and any 
specified memory locations Finally, it provides the means for system 
developers to modify the contents of memory, registers, status, or 
program counter, and to clear and set breakpoints. 
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refreshed. The chip then waits for a cycle when memory is 
not being accessed to do the refresh. There is a two-element 
queue in case a vacant memory cycle does not come before 
the following location must be refreshed. 

Keyboard Controller 

The circuit needed just to interface to the keyboard would 
have been quite small. Consequently, four programmable 
timers were added to this NMOS chip. The timers can count 
up to 27 hours with one-millisecond resolution. Each has a 
maskable interrupt. When the timer gets to its preset count, 
it interrupts, resets and counts again. 

The keyboard portion of the chip can scan an 8x 10-key 
keyboard plus three dedicated keys— shift, control, and 
caps lock. The key debounce time is mask programmable 
from 1.67 to 11.69 ms. After the key has been debounced. 
the chip generates an interrupt. An internal ROM is used to 
convert the key position to its ASCII* equivalent. 

Another feature of the chip is an output to a speaker for 
audio tones. Via firmware, a 1.2-kHz output tone can be 
obtained, or the frequency can be varied by periodically 
setting and clearing an internal flip-flop. This output can be 
accessed by the programmer using a command that can 
specify the frequency and duration of the tone. 

CRT Controller 

The NMOS CRT controller 1 interfaces the CPU to a 
127-mm diagonal CRT and its dedicated display memory. 
The memory is four 16Kx l dynamic RAMs. This is enough 
storage to hold a full display of graphics information 
(256x192 dots) plus four displays of alphanumeric infor- 
mation (4x32 characters/line x 16 lines/display). The chip 
must make sure that this memory is properly refreshed 
during the vertical retrace. The interface to the video drive 
circuitry consists of a vertical sync, pulse (60 Hz), a horizon- 
tal sync pulse (15.7 kHz), and a video line (4.9 MHz). An 
internal ROM provides a dot pattern translation of the 
stored ASCII characters. 

Printer Controller 

The NMOS printer controller 2 is the interface between 
the CPU and the thermal moving-head printer. To perform 
this function, the circuit must control a printhead drive 
motor, a paper advance motor and an eight-dot printhead. 
An internal 32-byte RAM allows firmware to buffer one 
print line of alphanumeric data at a time. An internal 
ROM translates an ASCII character into its appropriate 
dot pattern. 

Tape Controller and Read/Write Amplifier 

The NMOS tape controller and the bipolar read/write 
amplifier-' control the tape unit in the HP-85. The read/write 
amplifier reads and writes data on the tape through a two- 
track magnetic head. The tape is formatted using a 1:1.75 
delta distance code, in which 8-kHz flux reversals represent 
zeros and 4.6-kHz flux reversals represent ones. 

The tape controller encodes digital information into this 
format and sends it to the read/write amplifier. When read- 
ing the tape, the signals from the read/write amplifier are 
decoded by the controller. The tape drive motor direction 
and speed are also controlled by the chip. 

"American Standard Code (or information Inrercnange 



Buffer 

The NMOS buffer is the ninth chip in the system. With its 
help, the system bus can be expanded to include all of the 
plug-in I/O slots. It is capable of driving a 150-pF load. The 
delay between NMOS-level input and NMOS-level output 
is 50 ns. The signals that the buffer passes are eight bus 
lines, three control lines, and the interrupt and halt lines. 
The bus and control lines are designed to be bidirectional. 

The clock lines are not buffered since this would cause 
unwanted skewing. The clock generator is capable of driv- 
ing all internal as well as all external loads. 

Acknowledgments 

A great deal of credit must be given to the team that made 
this LSI chip set possible. Tim Williams designed the RAM 
controller and buffer. Donn Wahl designed the ROM. The 
keyboard controller was an effort from Jim Hutch ins and 
Jerry Erickson. Jerry also designed the CRT controller. The 
printer controller was the work of Clement Lo. Implementa- 
tion of the NMOS CPU was done by Jim Axtell. The car- 
tridge controller design is that of Doug Collins. The sense 
amplifier is the result of the combined efforts of Mike 
Moore. Doug Collins and Mike Barbour. Special thanks 
must be given to Don Hale and Howard Shishido, who did 
the majority of the chip mask designs. Howard Honig con- 
tributed to the project by writing most of the test programs 
for the LSI circuits. Mike Pan, Rock Davidson, and Tom 
Kraemer also contributed to the chip set design effort. Bob 
Tillman and Steve Larsen helped get the NMOS process and 
polysilicon resistor development up to production stan- 
dards. Payne Freret contributed some good ideas to the CPU 
design. Rick Bell did the design of the clock oscillator and 
generator. 

Reference 

1. J.F. Bausch. "A High-Quality CRT Display for a Portable Com- 
puter," Hewlett-Packard Journal, July 1980. 

2. C.C. Lo and R.W. Kelt, "A Compact Thermal Printer Designed 
for Integration into a Personal Computer." Hewlett-Packard Jour- 
nal. July 1980. 

3. D.J. Collins and B.G. Spreadbury, "A Compact Tape Transport 
Subassembly Designed for Reliability and Low Cost," Hewlett- 
Packard lournal, |uly 1980. 



Todd R. Lynch 

Todd Lynch is a native of Rochester, 
New York. He pined HP Laboratories in 
1 972 after working for a year on process 
control logic designs. Since coming to 
HP he has worked on computer ar- 
chitectures, and since transferring to 
what is now HP's Corvallis Division. 
, Todd designed the architecture for the 
HP-85 CPU He was the project man- 
ager for the integrated circuits in the 
HP-85 and is now project manager for 
high-end mainframes He received a 
BS degree from Grove City College. 
^ Pennsylvania in 1970 and an MSEE 

^Mk I from the University of Wisconsin in 

1971 Todd is co-chairman for the Willamette Valley Junior Tennis 
Tournament and in addition to tennis enjoys woodworking, softball. 
and skiing. He lives in Albany, Oregon with his wife and new son 




22 HEWLETT-PACKARD JOURNAL AUGUST I960 



© Copr. 1949-1998 Hewlett-Packard Co. 



Handheld Calculator Evaluates Integrals 



The HP-34C is the first handheld calculator to have a key 
that performs numerical integration almost automatically. 
It may change your attitude towards what used to be 
regarded as a dreary tedious task. 

by William M. Kahan 



NUMERICAL INTEGRATION has been the subject 
of about two thousand books and learned papers, 
with a dozen or so "new" methods published every 
year. And yet the task In question has a simple geometrical 
interpretation seen in Fig. 1: given an expression f|u) and 
lower and upper limits y and x respectively, the value 

I =/ y X fluldu 

represents the area under the graph of f(u| for u between y 
and x. Why so much fuss? 

As I write this an electrical engineering colleague. Pro- 
fessor J. R. VVoodyard, enters my office and asks to have 

' J„ 1 u-l In u ' 

evaluated on my IIP-34C Calculator (Fig. 2). Let's do it. 

Step L Key into the calculator under, say. label A a pro- 
gram that accepts a value u in the display (X regis- 
ter) and displays instead the computed value of the 
integrand 

V'u'(u-l) - l/|ln u) 

Ptg, 3 shows an HP-34C program to do this. 
Step 2. Restore the calculator to RUN mode and set the dis- 
play to, say, FIX 5 to display five decimal digits after 
the point, which areas many digits of the integrand 
as my client says he cares to see. (More about this 
later.) 

Step :J. Key in the lower and upper limits of integration 
thus.O ENTER' 1, thereby putting (I into the Y reg- 
ister and 1 into X. 

Step 4. Press Jy A. wait 25 seconds until the display shows 
0.03662. then press xSy to see 0.00001. We 
have just calculated 

I, = 0.03662 ± 0.00001. 

That was easy — too easy. Woodyard smiles as if he knew 
something I don't know. Could the calculator be wrong? 
How does the calculator know the error lies within ±0.00001 ? 

Many other questions come to mind: 

■ Why is numerical integration impossible in general? 
Why do we persist in trying to do it anyway? 

i How do we do it? How well do we do it? 

■ How does the Jy key compare with other integration 
schemes? 



' What can go wrong and how do we avoid it? 
• What else have we learned? 

These questions and others are addressed in the following 
pages. 

Tolerance and Uncertainty 

Integrals can almost never be calculated precisely. How 
much error has to be tolerated? The /y key answers this 
question in a surprisingly convenient way. Rather than be 
told how accurately / = jy f(u)du should be calculated, the 
HP-34C asks to be told how many figures of f(u) matter. In 
effect, the user is asked to specify the width of a ribbon 
drawn around the graph of f(u). and to accept in place of 
/ an estimate of the area under some unspecified graph lying 
entirely within that ribbon. Of course, this estimate could 
vary by as much as the area of the ribbon, so the calculator 
estimates the area of the ribbon too. Then the user may 
conclude from Fig. 4 that 

I = (area under a graph drawn in the ribbon) 
± IVi area of the ribbon) 

The calculator puts the first area estimate in its X register 
and the second, the uncertainly, in the Y register. 

For example. f(u) might represent a physical effect whose 
magnitude can be determined only to within, say, ±0.005. 
Then the value calculated as f(u) is really f(u) ± Af(u) with 
an uncertainty Af|u) = 0,005. Consequently FIX 2. which 
tells the calculator to display no more than two decimal 
digits after the point, is used to tell the calculator that 
decimal digits beyond the second cannot matter. Therefore 
the calculator need not waste time estimating I ± Af = 
Jy(f(u) ± Af(u))du more accurately than to within an 
uncertainty A/ = |/2 Af(u)du| This uncertainty is esti- 
mated together with / * AJ. thereby giving the calculator's 
user a lair idea of the range of values within which 1 
must lie. 









<(u> 




Area 1 = 


I'' t(u) du 





u=y u=x 



Fig. 1. An integral interpreted as an area 
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Fig. 2. HP-34C Calculator has keys to solve any equation and 
to compute integrals 

The uncertainty Af(u) is specified by the user via the dis- 
play setting. For instance. SCI 5 displays six significant 
decimal digits, implying that the seventh doesn't matter. 
The HP-34C allows the user's f-program to change the dis- 
play setting, thereby providing for uncertainties Af(u) 
that vary with u in diverse ways. But users usually leave 
the display set to SCI 4 or FIX 4 without much further 
thought. 

By asking the user to specify Af(u) instead of AJ the 
HP-34C helps avoid a common mistake — wishful think- 
ing. Other integration procedures, which conventionally 
expect the user to specify how tiny Ai should be. blithely 
produce estimates of / purporting to be as accurate as the 
user wishes even when the error Af(u) is far too big to justify 
such claims to accuracy. The HP-34C does not prevent us 
from declaring that f|u) is far more accurate that it really is. 
but our attention is directed to the right question and not dis- 
tracted by questions we cannot answer. Whether we spec- 
ify Af after a careful error analysis or just offer a guess, we 
get estimates J ±AI that we can interpret more intelligently 
than if we got only I with no idea of its accuracy or inaccuracy. 

A Survey of Integration Schemes 

Students are taught the fundamental theorem of calculus: 

J=f f(u|du = F(x)-F(y) provided F(u) = f(u) 
Jv du 

This means that one could calculate I if one could dis- 
cover somehow an expression F(u) whose derivative is 
the given expression f(u). Students are taught integration 
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as a process, applied to expressions, that starts with f and 
ends with F. But in professional practice that process 
hardly ever succeeds. A compact expression F(u) is al- 
most always difficult or impossible to construct from any 
given f(u). For instance, neither 

/* T exp(-u-/2)du»V27T nor /J exp(-u + x In u)du 

possesses a closed form, that is. an expression involving only 
finitely many elementary operations ( + , -, x, In. exp. tan. 
arctan. ...) upon the variable x. Nevertheless, both integrals 
can be approximated arbitrarily accurately by aptly chosen 
formulas. So often do statisticians and engineers need 
values of those integrals that formulas for them, accurate 
to ten significant decimal digits, ran now be calculated in 
a few seconds by pressing a key on certain handheld cal- 
culators. (Press Q on the HP-32E to get the first integral, 
the cumulative normal distribution: press x! on the HP-34C 
to get the second integral, the gamma function F(l+x), 
whether x be an integer or not,) 

Almost every rare integrand f(u) whose indefinite in- 
tegral F(x] = J* f(u)du + c is expressible in a compact or 
closed form can be recognized by a computer program that 
accepts the string of characters that defines f and spews out 
another string that represents F. (Such a program is part of 
the MACSYMA system, developed at MIT. that runs on a 
few large computers — two million bytes of memory — at 
several universities and research labs.) Perhaps the terms 
"compact" and "closed form" should not be attached to the 
expression F(x). since usually, except for problems as- 
signed to students by considerate teachers, the integral 




Begin with u in the X register 
Save u in the slack 

. . . \ u 
Recall u 

. . . \ u)(u-1) 
Recall u again 

. . . 1/ln(u) 



Display \ u/(u-1) - 1 In u 

Fig. 3. This program makes the HP-34C calculate the inte- 
grand Vul(u-1) - 1lln u when the argument u is in the X 
register and key A is pressed Labels B.O, 1 ,2. or 3 would have 
served as well as A 
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11 = %5 Area ot Ribbon = f" Al(u) du 




f(u) ? 



u = y 



Fig. 4. The graph o' an uncertain integrand t(u) -ArYu) can run 
anywhere in the ribbon bounded by the dashed lines The area 
under such a graph. /±A'. 'S uncertain by ±A'. which is one- 
half the area ol the nbbon The HP-34C displays its estimate of 
I'M m its K. register and holds an estimate ot M in its 1 register 

F(x) far exceeds the integrand f(u| in length and complex- 
ity- Shown in Fig. 5 are two compact forms and one closed 
form for F(x| when flul^l'll + u* 4 ). The extent to which F(x| 
is here more complicated than f(u) is atypically modest out 
of consideration for the typesetter. The formulas in Fig. 5 
will remind many readers of hours spent on calculus prob- 
lems, but they do not provide economical ways to calculate 
F(x) for any but very big or very tiny values of x. When I use 
the HP-34C's /y key to calculate F(l) = /□ du/fl+u 64 ) = 
0.989367 ± 0.000004 the answer appears in 200 seconds 
including 20 seconds taken to enter the f-program plus ISO 
seconds for a result (in SCI 5). Calculating F(l| from any 
formula in Fig. 5 takes at least about ten times longer, not 
including the time taken to deduce the formula. Fngineers 
and scientists have long been aware of the shortcomings of 
integration in closed form and have turned to other 
methods. 

Perhaps the crudest way to evaluate /J f|u)du is to plot 
the graph of l(u], like Fig, 1. on uniformly squared paper 
and then count the squares that lie inside the desired 
area. This method gives numerical integration its other 
name numerical quadrature. Another way, suitable for 
chemists, is to plot the graph on paper of uniform density, 
cut out the area in question, and weigh it. Engineers used 
to measure plotted areas by means of integrating engines 
Called planimeters. These range from inexpensive hatchet 
planimeters of low accuracy to Swiss-made museum 
pieces costing hundreds of dollars and capable of three 
significant decimals. (For more details see reference 1). 
Nowadays we reckon that the computer will drive 
the graph plotter SO it might as well integrate too, 

Today's numerical integration techniques are best ex- 
plained in terms of averages like 



A 



(x-y) =/ y x f(u)du/(x-y) 



which is called "the uniformly weighted average of f(u| 
over the interval between x and y." Another kind of average. 

it n 

A = 2 w,f( u,| where Wj>0 and 5)wj— 1, 
ii ii 

is .i finite weighted average ol n samples l'(u,|. f(u a ) I|u n |. 



Provided the sample arguments u,. u 2 u n . called nodes. 

all lie between x and y the sample average A will approxi- 
mate, perhaps poorly, the desired average A. and hence 
provide I = (x— y)A as an approximation to 1 = |x-y)A 
Statisticians might be tempted to sprinkle the nodes Uj ran- 
domly between x and y — that is what Monte Carlo methods 
do, But randomness is a poor substitute for skill because 
the error A - A tends to diminish like I V n as the number 
n of random samples is increased, whereas uniformly 
spaced and weighted samples provide an error that dimin- 
ishes like 1 n-. Other more artful methods do even better 
Different numerical integration methods differ princi- 
pally in the ways they choose their weights Wi and nodes tU 
but almost all have the following characteristics in com- 
mon. Each average A is associated with a partition of the 
range of integration into panels as shown in Fig. 6. Each 
panel contains one node u ( whose respective weight is 

Wj = (width of panel j)/|width of range of integration). 

The formula given above for A amounts to approximating 
the area in each panel under the graph of f(u) by the area of 
a rectangle as wide as the panel and as high as the sample 
f(Ui). The simplest method is the midpoint rule, whose 
nodes all lie in the middles of panels all of the same width. 
Other methods, like the trapezoidal rule and Simpson's 
rule, vary the panel widths (weights) and nodes in ways 
designed to exploit various presumed properties of the 
integrand f(u) for higher accuracy. Which method is best? 
If this question had a simple answer there would not be so 
many methods nor would we need texts like "Methods of 
Numerical Integration" by P.). Davis and P. Kabinowitz. 2 
which contains 16 FORTRAN programs and three bib- 
liographies with well over 10(10 citations. 

For example, consider Gaussian quadrature. This method 
is widely regarded as "best" in the sense that it very often 
requires fewer samples than most other methods to achieve 
an average A that approximates the desired A to within 
some preassigned tolerance. But the weights and nodes of 
Gaussian quadrature lake quite a while to calculate. Pro- 
grams to do so. and the resulting tables of weights and nodes 
for various sample counts n. have been published. 1 Had 
we chosen Caussiau quadrature for the Jy key we would 

F(x) = J f(u) du where t(u) - 1'(1 • u M ) 



F(x) = x£ (-x 6 V/(64k + 1) i! 
k o 



5 esc (£) sign Ix) - x V , x M )* / (64k - 1) it x> * 1 



^ V / Sln «. arctan ( . V,cos«, In 



1 t 



2 C0S«, 



-V 



■ (jjpj csc (gj) sign (x) il x ? 1) where", (k-Vi)ff 



132 



Fig. 5. Formal integration transforms many a simple expres 
sion l(u) into messy tormulas F(*) ot limited numerical utility 
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Fig. 6. The integral, regarded as an area, is here divided into 
lour panels each ol whose areas is approximated by the area 
of a rectangle as wide as the panel and as high as a sample. 

have had to store at least as many nodes and weights as we 
could expect to need for difficult integrals, amounting to 
at least several hundred 1 3-digit numbers, in read-only 
memory. But that would have left no space in the HP-34C 
for anything else, so a different method had to be found. 

The Jy key could not use a method that generates just 
one average A because that gives no indication of how 
accurately it approximates A. Instead we looked only at 
methods that sample repeatedly and with increasing 
sample counts nj < n 2 < n 3 < ... to produce a sequence of 

increasingly accurate averages A,, A 2 , A 3 Provided 

that sequence converges to A so fast that each | A^-i - A | 
is considerably smaller than its predecessor, the error 
|A^— A\ can be approximated accurately enough by 
| A^-Aj-t, |. and the last average A^.j can be accepted in 
lieu of A as soon as |A| t -A| c ^ 1 | is tolerably small. 

How small is "tolerably small"? That depends upon the 
area of the ribbon discussed above under "Tolerance and 
Uncertainty." Since the integral / Jy f(u)du inherits an 
uncertainty AJ = | Jy Af(u)du| from the uncertainty Af(u] 
in the integrand, so does A = //(x-y) inherit an uncertainty 
1A = AJ/|x-y|. which may be approximated by 

n 

AA = Y w,Af( U) ) 

in the same way as A is approximated by A. Indeed. A and 
AA can be computed together since they use identical 
weights and nodes. And so the sequence A t , A 2 . A 3 . ... is 
accompanied by a sequence of respective uncertainty esti- 
mates AA,, AA 2 . AA 3 Now "tolerably small" can be 

defined to mean "rather smaller than AA^,." 

The foregoing argument provides an excuse for accept- 
ing A^+i in lieu of A whenever two consecutive estimates 
A k and Afc+j agree to within AA^+j. but it provides no 
defense against the possibility that convergence is not so 
fast, in which case A k and A k ^, might agree by accident 
and yet be both quite different from A. The Jy key waits for 
three consecutive estimates A k , A k+] , and A k + 2 to agree 
within AA kt2 . Only the most conservative integration 
schemes wait that long. While this conservatism strongly 
attenuates the risk of accidental premature acceptance of 
an estimate, the risk that three consecutive estimates might 
agree within the tolerance and yet be quite wrong cannot 
be eliminated. Later, under "How to Deceive Everv Nu- 



merical Integration Procedure." some such risk will be 
proved unavoidable, but the risk now is so small that 
further attenuation is not worth its cost. 

The combination of ignorance with conservatism is 
surprisingly costly. Had we known in advance that A^ 
would be accurate enough we would have calculated none 
of the other averages. Instead, waiting for three consecu- 
tive averages to agree could easily cost some methods al- 
most 6.2. r ) times as many samples as if only A k had to be 
calculated, and more than that if the sample counts n t , 
n 2 . n 3 , ... are not chosen optimally. Fbl the/y key we chose 
n± = 2 —1 and we chose a method whose successive aver- 
ages each share almost half of the previous average's 
samples, thereby preventing the cost of ignorance from 
much exceeding a factor of 4. 

Memory limitations precluded the use of another family 
of methods known as adaptive quadrature. These methods 
attempt to distribute nodes more densely where the inte- 
grand f(u) appears to fluctuate rapidly, less densely else- 
where where f(u) appears to be nearly constant or relatively 
negligible. They succeed often enough that the best 
general-purpose integrators on large computers are adap- 
tive programs like Carl de Boor's CADRE; this and others 
are described in reference 2. Alas, adaptive programs con- 
sume rather more memory for scratch space than the 
twenty registers available in the HP-34C, 

What Method Underlies the /* Key? 

The IIP-34C uses a Romberg method: for details consult 
reference 2. Several refinements were found necessary. In- 
stead of uniformly spaced nodes, which can induce a kind 
of resonance or aliasing that produces misleading results 
when the integrand is periodic, the Jy key's nodes are 
spaced nonuniformly. Their spacing can be explained by 
substituting, say. 



into 



/ = 



f ff|«-|« 3 ) |(t-tf')sfo 



and distributing nodes uniformly in the second integral. 
Besides suppressing resonance, the substitution confers 
two more benefits. One is that no sample need be drawn 
from either end of the interval of integration, except when 
the interval is so narrow that no other possibilities are 
available, and consequently an integral like 



J 



sin u 
u 



won't hang up on division by zero at an endpoint. Second. 
I = Jy f(u) du can be calc ulated effi ciently when 
f(u) = g(u)V|x-u| or g(u)\ (x-ul(u-y) where g(u) is 
everywhere a smooth function, without any of the expedients 
that would otherwise be required to cope with the infinite 
values taken by the derivative f'(u) at u = x or u = y. Such 
integrals are encountered often during calculations of areas 
enclosed by smooth closed curves. For example, the area of 
a circle of radius 1 is 
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/„ \ u(4-u)du = 3.14159 ± 8.8x 10 

which consumes only 60 seconds when evaluated in SCI 5 
and only 110 seconds to get 3.141 592654 ±1. 4 x 10 ~ 9 
in SCI 9 

Another refinement is the use of extended precision. 
13 significant decimal digits, to accumulate the sums that 
define A^. thereby allowing thousands of samples to be 
accumulated, if necessary, without losing to roundoff any 
more information than is lost within the user's own f-pro- 
gram. The last example's 10 significant decimal digits of 
rr could not have been achieved without such a refinement. 

How Does the |y Key Compare with Other Integrators? 

What most distinguishes the HP-34C's Jy key from all 
other schemes is its ease of use. No step-size parameters, no 
plethora of error tolerances, no warning indicators that 
"can usually be ignored." Only the minimum informa- 
tion needed to specify Jy (f(u)~Af(u))du has to be supplied. 
And because the Jy key is effective over so wide a range of 
integrals it ranks among the most reliable procedures avail- 
able anywhere. Usually it is far faster than simpler proce- 
dures like the trapezoidal rule or Simpson's rule com- 
monly used previously on calculators. For integrands de- 
fined by programs that fit comfortably into a mid-sized 
handheld calculator that can hold at most 210 program 
steps, the/y key is comparable in speed (count the number 
of samples) with the integrators available on large com- 
puters. For much more complicated integrands the best 
adaptive integrators on large computers are appreciably 
faster. 

One of the HP-34(7s most important components is its 



Owner's Handbook. It is for most owners the first guide to 
the foothills of an awesome range of new possibilities. Two 
chapters are devoted to Jy. The first is introductory, and 
allows the user to evaluate simple integrals effortlessly 
and confidently. The second chapter is a longer explana- 
tion of the power and the pitfalls, concerned mainly with 
numerical integration generally rather than with the 
HP-34C in particular. This chapter had to be included be- 
cause its explanations and practical advice are not yet to be 
found in any text likely to be consulted by an owner, nor 
are they supplied by the instructions that accompany other 
integrators on other computers or calculators. This second 
chapter is part of the educational burden that must be 
borne by innovators and pioneers. The Owner's Handbook 
provides no formulas for the nodes and weights used by 
the HP-34C because they are not needed to understand 
how the Jy key works; instead they can be deduced from 
information in this article. 

Every numerical integrator like Jy. which executes a 
user-supplied program to get the integrand's value f(u|. 
imposes constraints upon that program. Some constraints, 
like requiring f to have a smooth graph on the interval of 
integration, are practically unavoidable. Others are nui- 
sances like 

■ Begin the f-program with a special label, say A'. 

■ Do not use certain memory registers, say #0 - #5. 

■ Do not use certain operations, say = and CLR. 

The /y key is encumbered with no such nuisances. The 
f-program may begin with any of several labels, so several 
different integrals can be calculated during one long com- 
putation. The f-program may use memory registers freely 
and may use any operation key except Jy itself. One of 
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Fig. 7. A program to evaluate 
I = Jo udu vtu) where u = v(u) 
satisfies v-u +in (1 *uve"]=0. 
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those keys is the HP-34("s powerful SOLVE key. 4 Conse- 
quently this calculator is currently the only one that can 
evaluate conveniently integrals of implicit functions. 
For example, let v = v(u) be the root of the equation 

v — u + ln(l + uue") ■ 0. 

Then 

j' t u du/u(u) = 1.81300 - 0.000005 

results from a program rather shorter than on any previous 
calculator; it is exhibited in Fig. 7. 

Furthermore, Jy may be invoked, like any other function, 
from within a program, thereby permitting the HP-34C to 
SOLVE equations involving integrals. For example, 
solving 

H cos (x sin «|d0 = 0 

for x = 2.405... takes a short program contained in the 
Owner's Handbook, and exhibits the first zero of the Bessel 
function l 0 (x). 

How to Deceive Every Numerical Integration Procedure 

Such a procedure must be a computer program — call it 
P — that accepts as data two numerical values x and y and 
a program that calculates f(u) for any given value u, and 
from that data P must estimate J = J*. f(u) du. The integra- 
tion procedure P is not allowed to read and understand 
the f-program but merely to execute it finilely often, 
as often as P likes, with any arguments u that P chooses. 
What follows is a scheme to deceive P. 

First ask P to estimate / for any two different values x and 

y and forf(u)s0. Record the distinct arguments . u 2 u n 

at which P evaluates f(u). Presumably when P finds that 
f(u,l = f(u 2 ) = ... = f(u n ) = 0 it will decide that I = 0 and 
say so. Next give P a new task with the same limits x and y 
as before but with a different integrand 

ffU] - |(u-u l )-|u-u,)-...-|u-u n |) i . 

Once again P will calculate f(ui). f(u 2 | and finding no 

difference between the new f and the old. P will repeat 
exactly what it did before. But the new integral / is quite 
different from the old. so P must be deceived. 

The HP-.'i4C's /y key can be hoodwinked that way. Try 
to evaluate/ J,' 28 f(u)du using first f(u) = 0 programmed 
in a way that pauses (use the PSE key) to display its argu- 
ment u. The calculator will display each sample argument 
it uses, namely 0. ±88. ±47 and ±117. Next program 

f(u| = |u(u-8B)|u + 88)(u-47)[u+-47)(u- 117)(u + 117)) 2 

and evaluate f-He flujdu again. The calculator will say 
that both integrals are 0. but the second polynomial's 
integral is really 1 .310269 x 10" 8 . That polynomial's graph, 
shown in Fig. 8, has the sharp spikes that characterize 
integrands troublesome for every numerical integration 
procedure. To calculate the integral correctly, reevaluate 
it as 2 J 0 128 f(u)du. thereby doubling the spikes' width com- 
pared with the range of integration. 



The threat of deceit impales the designer of a numerical 
integrator upon the horns of a dilemma. We all want our inte- 
grators to work fast, especially when the integrand f(u) 
is very smooth and simple like f|u) = 3u-4. But if the inte- 
grator is too fast it must be easy to deceive; fast integration 
means tew samples f(u } ). implying wide gaps between some 
samples, which leave room for deceitful misbehavior. 
Figs. 9a-9e illustrate the dilemma with two estimates of 
Jy f(u)du. The first estimate is based upon the three sam- 
ples drawn at the white dots, the second upon seven sam- 
ples including those three white plus four more black dots. 
Fig. 9a shows why all sufficiently smooth graphs f(u| that 
agree at all seven samples have nearly the same integrals, 
but Fig. 9b shows how two integrands could provide the 
same samples and yet very different integrals. The coinci- 
dence in Fig. 9b is unlikely; successive estimates based 
upon increasingly dense sampling normally would reveal 
the difference as in Fig. 9c. However, situations like those 
illustrated in Figs. 9d and 9e are very likely to deceive. 

Textbooks tell us how to avoid being deceived: avoid 
integrands f(u) among whose first several derivatives are 
some that take wildly different values at different places 
in the range of integration. Or avoid integrands f(u| that 
take wildly different values when evaluated at complex 
arguments in some neighborhood of the range of integra- 
tion. And if wild integrands cannot be avoided they must 
be tamed. We shall rejoin this train of thought later. 

Improper and Nearly Improper Integrals 

An improper integral is one that involves in at least 
one of the following ways: 

i One or both limits of integration are ±f-, e.g., 

1.8 

X -rr 

10" 



10"-- 

II II 
II II 

!(u) = (u(u 2 -47*)(u J -8B J )(u J -117 J )) J 




128 -88 0 47 117 128 



Fig. 8. The polynomial Ifu) was devised to deceive the 
HP-34C into miscalculating its integral as 0 instead ol 
1.31 x JO 28 This spiky graph is typical ol integrands that can 
battle any numerical integrator 73% ot the area under the 
graph lies under t wo spikes whose widths span less than 4 % ot 
the area ot integration 
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J_,exp f-u J du m Virf2. 



- 




Fig. 9. Few samples (open circles) mean last integration but a 
large possibility ot error More samples (solid dots plus open 
circles) usually mean more accuracy, but not always , as in (b), 
(d), and (e) (a) Which is the graph ol Hup No matter, both 
have almost the same integral, (b) Which is the graph ol f(up 
They have very ditlerent integrals, (c) Here two graphs that 
coincide on the first samples O are distinguished by a signifi- 
cantly ditlerent outcome after second samples • are drawn (d) 
It the graph of t(u) has a few sharp and narrow spikes, they will 
probably be overlooked during the estimation ot the integral 
based on finitely many samples, (e) II the graph ot t(u) has a 
step that was not made known during the estimation of the 
integral, then the estimate may be mistaken 



* The integrand tends to ±* someplace in the range of 
integration, e.g.. 

/„' ln(u)du = 1. 

■ The integrand oscillates infinitely rapidly somewhere in 

the range of integration, e.g.. /,} cos (In u)du = 
Improper integrals are obviously troublesome. Equally 
troublesome, and therefore entitled to be called nearly 
improper, are integrals afflicted with the following malady : 
= The integrand or its first derivative changes wildly with- 
in a relatively narrow subinterval of the range of integra- 
tion, or oscillates frequently across that range. 
This affliction can be diagnosed in many different ways. 
Sometimes a small change in an endpoint can render the 
integral improper, as in 

JTq'ooo, ln(u)du =- 0.99898 ... -*£ ln(u)du = 1, 

Sometimes a small alteration of the integrand can render 
the integral improper, as in 

j\dxl[x 2 + Hr'") = 314157.2654... -»/' dxftc* = x. 

Sometimes the value of the integral is nearly independent 
of relatively huge variations in one or both of the end- 
points, as is /„ exp (-u 2 ) du = Vw/2 for all x > 10. Regard- 
less of the cause or diagnosis, nearly improper integrals 
are the bane of numerical integration programs, as we 
have seen. 

During the HP-34C's design a suspicion arose that most 
integrals encountered in practice might be improper or 
nearly so. Precautions were taken. Now that experience 
has confirmed the suspicion, we are grateful for those pre- 
cautions. They were: 

1. Avoid sampling the integrand at the ends of the range of 
integration. 

2. By precept and example in the Owner's Handbook, 
warn users against wild integrands, suggest how to recog- 
nize them, and illustrate how to tame them. 

The second precaution ignited controversy. Against it on 
one side stood fears that its warnings were excessive and 
might induce paranoia among potential customers. Who 
would buy a calculator that he thinks gets wrong answers? 
Actually wrong answers were very rare, thanks in part to 
the first precaution, and many attempts to vindicate dire 
predictions about mischievous improper and nearly im- 
proper integrals were thwarted by unexpectedly correct 
answers like 

j' t In (u) du = 0.9998 ± 0.00021 
in 2 minutes at SCI 3. Or 

■> 

J expl-u") du = 0.886227 ± 0.0000008 

in 4 minutes at SCI 5. II the wages of sin be death. O Death, 
where is thy sting? 

On the other side stood a number of embarrassing ex- 
amples like 
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J 4 0 °° exp(-u 2 ) du 

miscalculated as 0.0 ± 0.0000000005 in 14 seconds. 
Another, had we known it then, would have been Wood- 
yard's example at the beginning of this article; the correct 
answer 

f ( — " r 2 - ) du = 0 03649 ± 0.00000007 
Jo " -1 In u * 

in 23 minutes at FIX 7 differs from FIX 5's wrong answer 
0.03662 in the worst way; the error is too small to be ob- 
vious and too large to ignore. Adding to the confusion 
were examples like 



0.12-r 



Alx) = 



J" V -2 In cos(u-] du/u- = 1 - x^O + x 8 '480 



for which computation in SCI 4 produced ridiculous values 
like 4(0.1) = 0.95742 ± 0.00005, 4(0.01) = 0.58401 ± 
0.00003. and 4(0.001) = 0. all impossibly smaller than 1. 
This example ap pears to conde mn the /y key until the 
integrand f(u) = v' -2 In cos(u 2 )/u 2 is watched for small 
arguments u and seen to lose most of its figures to round- 
off, losing all of them for I u I S0.003. despite an absence of 
subtractions that could be blamed for cancellation. Then 
the example appears to condemn the whole calculator. Who 
wants responsibility for a calculator that gets wrong 
answers? 

Don't panic! The answers are wrong but the calculator 
is right. 

How to Tame a Wild Integral 

Forewarned is forearmed. Every experienced calculator 
user expects to encounter pathological examples like some 
of those above, and expects to cope with them. The ques- 
tion is not "whether" but "when"? And that is when atten- 
tion to detail by the calculator's designers is rewarded by 
the user's freedom from petty distractions that can only 
complicate a task already complicated enough. But like 
the dog that did not bark,* the absence of distracting de- 
tails may fail to be appreciated. That is why the examples 
explained below have been chosen — to illustrate the ad- 
vantages of liberated thought. Work them on your calcula- 
tor as you read them; don't skim them like a novel. Then 
you may come to think of your calculator the way I think 
of mine, as a trusted friend who stays with me when I 
need help. 

The integral 4 (x) above contains an integrand f(u) = 
V —2 In cos(u 2 )/u 2 that loses its figures when u becomes 
tiny. The problem is caused by rounding cos(u 2 ) to 1, 
which loses sight of how small u 2 must have been. The 
solution compensates for roundoff by calculating f(u) 
as follows: 

Let y = cos u" rounded. 
If y = 1 then let f(u) = 1 
else let f(u) = V 



-2 In y'cos 1 y. 



The test for y = 1 adds four steps to the f-program and. 
provided In and cos -1 are implemented as accurately as 
on all recent HP calculators, the problem goes away. 

•See the last tew paragraphs ot the Snertock Holmes story Silver Bia^e" by Conao Doyle 
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Fig. 10. Substituting w 2 lor u turns the wild graph (a) into the 
easy one (b). But do not replace (w + 1)) by 

because roundofl errors introduce a spike, as shown in (c) 

Woodyard's example 1, has an integrand f(u) whose de- 
rative f'(u| — *« as u — »0. The graph of f(u) shown in Fig. 1 0a 
looks like a lovers' leap. Stretching the u-axis near u=0 by 
substituting u = w 2 turns the precipice into the hummock 
shown in Fig. 10b and transforms the integral into an easy 
calculation: 



n — ^ 

Jo 1 (w-l)(w + l) 



In w 



The HP-34C computes this as 0.03649 ± 0.000005 in 100 
seconds at FIX 5 or 0.0364900 ± 0.00000008 in 200 seconds 
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at FIX 7. Do not replace (w-l)(w+l) by (w 2 -l) because 
the latter loses to roundoff half of its significant digits as 
w — » 1 and introduces a gratuitous spike into the integrand's 
graph shown in Fig. 10c. which was plotted on an HP-85. Do 
not worn' about w = 0 or w = 1 because they don't happen, 
but do worry that as w — » 1 the integrand approaches the un- 
reliable expression » - * = 0. This means that FIX 7 dis- 
plays about as many digits as could possibly be correct for 
all w < 0.999. beyond which the JJ key draws few if any 
samples because it converges so fast. 

The graphs of exp(-u 2 ) over 0 s u s 300 and of 
l'(u 2 + 10 ~ 10 ) over -1 « u « 1 both resemble huddled mice 
with very long tails stretched out hundreds or thousands of 
times as long as iheir bodies. Plotting the graphs on a page 
of normal width is futile because the bodies get squashed 
into vertical whiskers. 

Most people who integrate such functions numerically 
cut off the tails. Thin tails can be cut almost indiscrimi- 
nately without much degrading the accuracy or the speed of 
integration. Such is the case for J Q exp (-u 2 )du, which J y 
evaluates in less than. say. 4 minutes at SCI 5 provided that 
x. if bigger than 4 or 10. is cut back to something between 4 
and 10. But /_* du/(u 2 + 10 ~ 10 ) has too thick a tail to cut 
without losing accuracy or patience when x is large. That 
is why Procrustean methods are not recommended. Better 
to shrink the tail via an artful substitution like u = K + /utan v 
where A lies within the body of the mouse and /x is roughly 
that body's width. Doing so with a = 0 and y. = 1 changes 
/ ft exp|-u 2 ) du into 

fwsun % >) . 

J 0 exp| - tan"v) (1 + lan _ v) dv 

which jy evaluates in three minutes ;il SCI 5 even when 
x is as big as 10 10 . Don't worry about tan w/2 because it can't 
happen on a well-designed calculator. 

J^du/lu 2 + in '"I 

benefits miraculously from the foregoing substitution 
when A = (J and n = 10 5 . but values near those do almost 
as well. 

Another technique might be called "subdivide and con- 
quer." It subdivides the range of integration into subinter- 
vals upon each of which the integrand f(u) is tame, al- 
though f{u) may l ook wild on the range as a whole, For 
example, f(u) = V u z + 10 ~ 10 has a V-shaped graph prac- 
tically the same as that of lul. Evaluating /.'' 3 f(u) du 
accurately takes a long time if done with one press of 
jy, but subdividing the integral into 

j" :) f(uidu +/,;f(uid U 

takes two presses of and one of £+ but much less time. 

Subdivide and conquer works best when combined with 
apt substitutions. For example, if the formulas in Fig. 5 
were unavailable how would Ff» = du/(l + u 64 ) be 
calculated? 

F(«o) = /J du/(l + u M ) + /* du/|l + u M ) ... subdivided 
= du/|l + u M ) + /„' Y,"-dwHw M + 1) ... u = l/w 



m Jj (1 + u 6; | du/|l + U M ) 

= ' - Jo" (u"-u M ) dum - u M ) 

= 1 -Vfc J 0 ! (1-V*| i> 5M dv/(l + v") 

= 1 000401708155 + 1.2X10" 1 * 



.. merged via w=u 

... some algebra 

... u=v l * to shrink 
a tail 



in 10 minutes at SCI 8. Thus we have calculated F(x) = 
(7r/64)csc(jr/64) to 13 significant decimals on a ten-sig- 
nificant-decimal calculator. 

Oscillatory integrals like /q cos(ln u) du sometimes suc- 
cumb to stretching substitutions like u = v 2 that damp the 
oscillations, but generally oscillatory integrals cannot be 
calculated accurately and quickly without sophisticated 
tricks beyond the scope of an article like this. A simple 
trick worth trying when the period of oscillation is known 
in advance is called folding, though it is really another 
instance of subdivide and conquer. Here is a didactic 
example. 



r 



^sin u_^ __ j u _ s( j|| runn i n g a f, er over three 
Vu+Vu+w hours at sci 5 



Vu + Vu-Hr 



du 



-2T 

,y n •* II 



Vu+n7T + Vi> + n7TT7r 



dv 



after being subdivided and with u = v + nTr. Exchanging 
/ and v produces 



l 

I] = I sin 2 v • 2. , 



dv. 



+ n7T + 7T 



At this point a program should be written to calculate 
the sum. but because the example is didactic the sum 
collapses to yield 



h = 



600 sin-v 



Vv+600tt 



dv = 21.10204 ± 0.00007 



in 5 minutes at SCI 5, 

Now for a final example drawn from life: 



V = 



du 



ii (a 2 +u) vV+u) (b 2 +u) (c 2 +u) 



fora = 100. b = 2,c = 1. 



This integral pertains to the electrostatic field about an 
ellipsoidal body with principal semiaxes a. b. c. 5 The 
ellipsoid is needle-shaped like an antenna or a probe. The 
classical approach transforms V into a standard form called 
an elliptic integral of the second kind and interpolates on 
two variables in published tables to get a numerical value. 
The following approach takes less time. 

First transform the improper integral (Jo) into a proper 
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one by substituting, say. u = (a 2 -c 2 ]/(l -u 2 ) - a 2 to get 
where 



V = A J V (l-u -)/(.,-+■ a) dv 



A = 2/1(8"-^) Va Z -b 2 ) = 2.00060018 x 10"" 

H = c'a = 0.01 

a = [tf-c'Vtf-b 1 ) = 3.001200480 X 10 1 
i 

Now, as always happens when a >> b > c. the integral 
is nearly improper because a and (i are both so nearly 0. We 
suppress this near impropriety by finding an integral in 
closed form that sufficiently resembles the troublesome 
part of V. One candidate is 

W=X / dv/vV+a = A In |v+vV+a) I ' 



Then 



= A In ((1+ Vl+al/l^ + vV 2 -*-")) 

= 8.40181880708X 10~ K 
V = W + A/' (V[\-v 2 )l[v 2 +a) - JIVv 2 +a) dv 



W/A 



-Ji (1 + Vl-u 2 ) Vv- + a 



) dv 



7.78867525X10" 6 ± 1.3x10 14 



after seven minutes at FIX 8. Don't worry about VI— v z 
as v * 1 because the figures lost to roundoff are not needed 
and its infinite derivative doesn't bother the HP-34C. 

Conclusion 

A powerful mathematical idea has been placed at the 
disposal of people who will invoke il with fair confi- 
dence by pressing a button marked Jy without having to 
understand any more about its internal workings than 
most motorists understand about automatic Iransmissions. 
Integrals that might previously have challenged the 
numerical expert and a big computer now merely amuse 
the scientist or engineer, and tomorrow they will be rou- 
tine. And now those engineering students who do attend 
classes in numerical analysis need no longer be expected 
to memorize the names nor the remainder terms of quadra- 



ture formulas but may instead be laught lo use integra- 
tion wisely. 
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